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Fellows of the Physical Society will realize that during the last few years 
the size of the Proceedings has significantly increased. The Papers Committee 
and Council have had to decide that this continued expansion must now be 
checked, both because of increases in the cost of paper, of printing and of 
the other services ancillary to publishing, and also because of the increasing 
difficulty of obtaining, for example, sufficiently large allocations of paper. 

In the matter of actual presentation of papers assistance can be given by 
the way in which diagrams and mathematics are prepared. Diagrams should 
be displayed in such a way as to take a minimum of space, e.g. a shallow, wide 
diagram is preferable to a long, narrow diagram, and individual single curves 
can often be superposed in one framework. In general no elementary 
mathematical steps should be included in a paper, and every effort should be 
made to reduce tables to a minimum, and to replace any series of values, 
wherever possible, by the mean and a standard deviation. ‘The co-operation 
‘of Fellows of the Society can materially assist in procuring these small but 
necessary economies. 

Papers are still sometimes written in an unduly verbose style. If the 
teferees of such a paper have made no objections and the paper is accepted, 
its undue length makes inroads into the publication space available. If the 
referees do object, the typing of referees’ reports and additional communications 
to authors and referees places an extra burden on the secretarial staff. If authors 
ensured that their papers were written in a concise style, considerable advantage 
would result: submitting a manuscript to a critical colleague can be a useful 
check. 

Recently an increasing tendency has been noticed for authors, both individual 
and in research establishments, to submit for publication papers which contain 
interim experimental results, or tentative conclusions from their work. It is 
felt that the value of such articles would be greatly increased if their publication 
were withheld until more final results and conclusions were established. 
Moreover, much valuable time on the part of the referees of the papers, and 
of the office staff, could be saved if authors would consider this point more 
seriously before submitting papers for publication. 


H. H. Hopkins, 
Honorary Papers Secretary. 
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ABSTRACT. Ammonia—oxygen, ethylene-oxygen and other diffusion flames were 
examined spectroscopically using the flat flame technique described in a previous paper. 
Each flame is shown to consist of a main reaction zone where the temperature is a maximum. 
The emission from this zone is thermal, and the component gases are in chemical 
equilibrium. On either side of the main reaction zone there is a region with a steep 
temperature gradient which constitutes a pre-heating zone. Fuels such as hydrocarbons 
or ammonia and oxidizers such as nitric oxide undergo chemical change in these 
pre-heating zones. On the fuel side the change is principally a thermal decomposition 
because the main reaction zone presents a barrier to oxygen penetration. ‘Thus carbon 
formation in hydrocarbon diffusion flames is essentially a pyrolytic process. The 
application and value of absorption spectroscopy are discussed and quantitative 
measurements of the concentrations of molecular oxygen and OH radicals are recorded 
which could not have been otherwise obtained. These results and a number of qualitative 
observations clearly show that diffusion flames with oxygen at atmospheric pressure have 
the structure predicted for a flame in which diffusion is the rate determining process. 
A number of differences are described between the characteristics of diffusion and 
pre-mixed gas or Bunsen flames. 


§1. INTRODUCTION 

NEW technique for the investigation of diffusion flames which was 

developed by the authors has been described in a previous paper 

(Wolthard and Parker 1949). The principal feature of the technique is 
the very full use of absorption spectroscopy which can be made with it. It has 
now been applied in a more quantitative way and from a detailed study of one 
particular problem, viz. the diffusion flame with oxygen, a number of new facts 
have been obtained which are considered in relation to the structure of the 
flame. ‘These facts help to show the essential differences between diffusion 
flames and those of the pre-mixed gas or Bunsen type. 


§2. DISCUSSION OF THE MECHANISM OF DIFFUSION FLAMES 

In order to obtain a mental picture of the physical processes in a diffusion 
flame let us suppose that we have two gases such as hydrogen and oxygen 
enclosed in a box and separated by a diaphragm. If the diaphragm can be 
removed and the gases lit instantaneously at their interface then the temperature 
will rise at this interface by the production of heat. After a very short interval 
of time the gases close to the interface will reach equilibrium and will have a 
definite composition of H,O, H,, H, O., O and OH, which will not change 
unless the temperature changes. ‘T'wo processes will now occur: 

(a) Heat will flow towards the unburnt hydrogen and oxygen, and will 
raise the temperature of these gases. The velocity with which this heat travels 


. 
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is governed by known laws of conductivity and is essentially an error function, 
_ Le. the heat travels rapidly at first on account of the steep temperature gradient 
one ear the source and then more slowly as the distance increases. Since the 
thermal conductivity changes markedly with temperature it is not feasible to 
onl culate this heat flow. 

(6) The unburnt hydrogen and oxygen will diffuse into the zone of burnt 
gases where their partial pressures are low. This diffusion process will be 


because the diffusion coefficients increase with increasing temperature. 

_ The entry of fresh gas into the zone of burnt gases would raise the partial 
ressures of the hydrogen and oxygen above the equilibrium amount. At the 
igh temperature involved however the velocity of the chemical reaction is so 
at that equilibrium is achieved at once and the heat released in the process 
places the heat lost by conductivity. 

_ Thus our conception of a diffusion flame is that of a mass of hot burnt gas 
formed at the centre. Diffusion, which is a relatively slow process, tends to change 
he composition of this gas but chemical reaction follows very rapidly and does 
allow any chemical potential (i.e. non-equilibrium condition) to develop. 
he mass of hot burnt gas grows thicker and the gradient of partial pressures 
ff the hydrogen and oxygen away from the centre becomes less steep. The 
overall velocity of combustion therefore slows down until the final flame 
femperature is reached and the partial pressure of each constituent of the burnt 
yas is uniform. 

Gases at the centre which are at 2000°c or higher showld begin to radiate 
especially if carbon monoxide, oxygen and OH molecules are present. Some 
ization will take place and the recombination of ions and electrons will cause 
ission of weak continuous light in the visible region of the spectrum. ‘The 
ht emitted by the flame comes from regions of sufficiently high temperature 
gas radiation and does not necessarily indicate the point of combustion. 
this respect diffusion flames contrast with pre-mixed gas flames where the 
iation of the inner cone comes directly from the reacting molecules. Thus the 
H, C, and CH radiations of the inner cone of a pre-mixed hydrocarbon—air 
flame are not those of molecules in thermal equilibrium (Gaydon and 
Wolfhard 1948). 

_ The conditions on ignition of the gases in the hypothetical case of the box 
brought about in effect in the flat diffusion flame mentioned previously. 
‘the burner (Wolfhard and Parker 1949) is supplied with hydrogen and oxygen 
carefully controlled rates, two parallel streams of these gases emanate with 
ual velocities and laminar flow from the two internal rectangular tubes of the 
When the flame is lit a series of events similar to those described for 
box occurs in i a vertical direction. Thus a stationary flame is obtained 


not the reaction velocity is the rate determining process in the flame. ‘This 
umption is believed to be valid for the following reason. 

_ A pre-mixed stoichiometric flame of ethylene and oxygen has a flame 
city of about 600cm/sec and a reaction zone of only (0-005 cm thickness. 
Thus the mean rate of conversion is about 4 mole of oxygen burnt per cm* 
I= 2 
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per second. The final flame temperature which may be calculated is only | 
reached after the gases have passed through the inner cone of the flame. 
It appears therefore that this very high reaction velocity is achieved in a part 
of the flame where the temperature is probably below the final temperature. 

In the diffusion flame combustion is rapid at first and then slows down until 
everything has been burnt but the reaction always takes place at a temperature 
close to the final temperature. At this temperature the activation energy of 
the reaction can be overcome and a high concentration of free radicals is 
presumably present to assist combustion. In spite of these advantages the 
reaction velocity is small compared with that of the pre-mixed gas flame. For 
example the rate of combustion in a diffusion flame of ethylene and oxygen was 
measured on our burner. With 50cm?/sec of each gas emanating from the 
two parallel tubes (each 5cmx1cm) the flame had a height of 7cm and a 
reaction zone about 1cm thick. The mean rate of conversion was therefore 
of the order 6 x 10-> mole of oxygen burnt per cm® per second. This rate is 
nearly five orders lower than that of the corresponding pre-mixed gas flame 
alhtough conditions of temperature and radical concentration should be’ 
favourable to the diffusion flame. It must be concluded therefore that the 
rate of combustion is not fixed by the reaction velocity but by the rate of diffusion 
of the reacting molecules. In such cases it is probably justifiable to regard the 
diffusion flame as a zone of hot gases in chemical equilibrium where diffusion 
processes create small concentration and temperature changes which are very 
quickly brought back into equilibrium by chemical reaction. 

This conception of the diffusion flame only holds in ideal cases, however, — 
such as those of the oxygen flames investigated here. In flames which have a 
lower temperature and a lower reaction velocity pre-mixing may occur in the 
reaction zone and some chemical potential may develop there. Diffusion will 
then cease to be the sole rate determining process. It is possible that air flames 
will come into this category but the question has not yet been investigated. 
Low pressure conditions may also change the relative rates of diffusion and 
chemical reaction. Thus if we imagine that the pressure of hydrogen and 
oxygen in our box is only 1 mm Hg instead of 760 mm Hg the reaction velocity 
assuming a second order reaction will be only (1/760)? but the diffusion rate 
will be 760 times greater. Mixing by diffusion will be almost instantaneous 
whereas reaction will take an appreciable time. We shall have in fact a pre-mixed 
gas flame. Experiments in support of this argument are mentioned later. 


|| 


§3. APPLICATION OF ABSORPTION SPECTROSCOPY TO 
DIFFUSION FLAMES 

Up to the present spectroscopic studies of flames supported on a round tube 
have been limited by the geometry of the flame. Thus the important reaction 
zone presents a three-dimensional problem in which it is impossible to locate 
the emitters exactly or to apply absorption spectroscopy satisfactorily. The 
emission spectrum which is obtained is only a measure of the concentration 
of the unexcited atoms or molecules when these are in thermal equilibrium and 
there is no self-absorption. Moreover the correlation of the strength of emission 
and concentration is a difficult matter, and can only be used in a few special cases 
to get an approximate value of the number of unexcited molecules. In the flat 
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diffusion flame the reaction zone at ordinary pressure is 5-10 mm thick and has 
an optical depth of 5 cm or more. ‘The thickness of the reaction zone may be seen 
in the photographs of the end-view of the ammonia and ethylene flames, figs. 1 
and 2 respectively. It is therefore practicable to get absorption as well as emission 
spectra of reacting molecules. This is an important advance since the strength 
of the absorption gives a direct measure of the number of molecules involved 
and thus the concentration, as well as the exact location, of some of the molecules 
may be obtained. 

As regards the general optical arrangement which is used for investigations 
with flat flames, it is essential to have quartz-fluorite achromatic lenses otherwise 
it is impossible to focus light of all wavelengths from the flame on to the slit of 
the spectrograph. The lens used between the flame and the spectrograph had 
a focal length of 40 cm and a diameter of 2-4 cm. It was stopped down to 0-8 cm 
in a direction parallel to the direction of gas flow. A very small aperture was 
employed in order that the different light zones of the 5 cm long flame would. 


Figure 1. Photograph of NH;—-O, Figure 2. Photograph of C,H,-O, 
flame; end view, actual size. flame; end view, actual size. 


not overlap appreciably in the image. By placing the lens 80 cm from the flame 
unit magnification was obtained and the overlapping in the image of the flame 
was not greater than +0-01 cm. 

In absorption experiments a spurious absorption effect was noticed which 
was caused by refraction. Light from the background light source had to pass 
through the flame and therefore through regions with steep temperature gradients. 
In doing so it was refracted (as in a Schlieren system) away from the lens and was 
partially or wholly missing from the image on the slit. ‘The impression of a 
continuous absorption was thus obtained from the photographic plate. This 
effect was troublesome at first until the reason was established. Fortunately, it 
only occurred outside the main flame in regions of lower temperature. Schliere 
develop in regions of rapid change in refractive index, i.e. where 1/7 changes 
rapidly in hot gases, and are therefore strongest at rather low values of the 


absolute temperature 7. 
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Interpretation of Absorption Spectra 


The spectrographs used gave strong and measurable absorption but did not 
truly resolve the lines of band spectra. The measured strength of such an 
absorption is a complicated function and depends mter ala on whether the line 
width consists mainly of Doppler broadening or of collision broadening. Van der 
Held (1931) has considered the problem and has concluded that for pure collision 
broadening the absorption strength is proportional to the number of absorbing 
centres. For pure Doppler broadening on the other hand, the strength of the 
absorption increases proportionally at first but quickly approaches a saturation 
value. Practically there is a combination of these two effects, and in the present 
investigation it was found that the absorption increased proportionally at first 
but that the rate of increase subsequently became smaller with increasing numbers 
of absorbing centres. 

These results indicate that the absorption coefficient is not proportional to 
the number of absorbing centres, i.e. Beer’s law does not hold. To overcome 
this difficulty we compared the strength of absorption in the flame under investi- 
gation with the absorption in another flame where conditions were known. 
A number of precautions were necessary which are best described in relation to 
a particular system. In the OH system we used the (0,0) band with its head at 
30644. The absorption is only slightly affected at different temperatures by the 
changing population in different vibrational levels because most of the molecules 
are in the lowest vibrational state at the temperatures encountered. Within this 
band some Q lines in front of the second head (Q 34, 43, 53, 64) are most sensitive. 
The strength of absorption for the rotational levels changes with temperature 
and the change in population may be calculated. ‘These changes inside the 
rotational population were small for the temperature range covered by our 
measurements, viz: 2200°c-2800°c, and could be disregarded because the 
overall accuracy of measurement of the OH and O, concentrations was estimated 
to be about 20-30%. 

The strength of the absorption could be calibrated in terms of the OH con- 
centration by the use of a pre-mixed gas flame of similar optical depth. Pre-mixed 
gas flames are in chemical equilibrium above the inner cone and a pre-mixed gas 
flame of ammonia—oxygen or ethylene-oxygen was burned on a thin burner of 
about 0-3 mm width which had a length of 5 cm, similar to the optical depth of 
the flat diffusion flame. The flame velocity of the pre-mixed gases is high and 
it is reasonable to assume that at a distance of 1-2 mm above the inner cone where 
the flame width is 8-10 mm the temperature reaches the theoretical flame 
temperature corresponding to the mixture strength. The theoretical flame 
temperatures and equilibrium concentrations of the gaseous products corres- 
ponding to the temperature are shown in figs. 7 and 8 plotted against the mixture 
strength. ‘The latter is expressed as the ratio n,:,. It has been assumed that 
there is no external heat loss and it was also ascertained that no undecomposed 
ammonia passed through the ammonia-oxygen flame even at very low ny: ny 
ratios. Similarly in the ethylene-oxygen flames no carbon formation occurred 
in the n):m, range covered by the graph of fig. 8. 

At low n,:n, ratios, i.e. in fuel-rich mixtures, the OH concentration rises 
steeply towards the stoichiometric ratio and falls again more gently as the n,: m, 
ratio increases above stoichiometric. Mixtures in the region of the stoichiometric 


and on the fuel-rich side were therefore used for getting a calibration curve of 


absorption strength versus the OH concentration. 'The absorption strength 
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was taken as the ratio of the intensity in the absorbing line to the intensity of the 
unabsorbed continuum (of the light source) beside the line. The intensity in 
each case was estimated from the density of the plate after the latter was calibrated 
in the usual way. A Hilger large Littrow-type spectrograph was used for the 
OH system in conjunction with the anode of the carbon are as background 
illumination. 

In comparing the strength of the OH absorption from the pre-mixed flame 
with that in the diffusion flame, it is necessary to have a similar gas composition 
(i.e. H,O, H,, O, etc.) in each case because these gases affect the collision 
broadening of the OH. This was done as far as possible, though small variations 
would presumably cause only small errors since collision broadening itself does 
not contribute much to the width of the line in very hot flames. 

Measurements of oxygen concentration were more restricted than those 
of the OH molecules. Molecular oxygen cannot be measured in its (0,0) 
Schumann—Runge band as this has a low transition probability and is outside 
the range of the quartz spectrograph. Hot oxygen absorbs, however, because 
higher vibrational levels are excited which permit a transition to the upper 
electronic state above 20004. Both pre-mixed and diffusion flames using 
oxygen show strong oxygen absorption in the region 2 200 A and in cases of excess 
oxygen this absorption may be followed up as high as 26004. Although this 
absorption has been observed before (Fiichtbauer and Holm 1925), a satisfactory 
analysis has not been achieved up to the present. The strength of the absorption 
clearly depends on the temperature of the flame and it is advisable to measure 
as far as possible in the far ultra-violet region to minimize the temperature effect. 
Information on collision broadening is not available. ‘The concentration of 
oxygen molecules in pre-mixed gas flames is only appreciable for lean mixtures 
and in the diffusion flame a high concentration occurs only on the oxygen side 
of the main reaction zone. The actual calibration of the strength of the absorption 
against the oxygen concentration was done in the manner described above using 
oxygen-rich pre-mixed gas flames. It was fortunate for the accuracy of the 
measurements that the temperature in the hottest part of the diffusion flame is 
close to the temperature of the stoichiometric pre-mixed gas flame and the 
strengths of the O, absorption seemed just about to match. Moving from the 
point of highest temperature outwards on the oxygen side the O, concentration 
increased while the temperature decreased. Pre-mixed gas flames had a similar 
trend as the mixture strength was weakened off. On the outside of the diffusion 
flame where the oygen was cold there was no absorption, and no comparison could 
be made. The main interest of this work, however, centred round the high tem- 
perature region of the diffusion flame. A check on the method was afforded by 
using two oxygen lines of different wavelength (corresponding to different 
vibrational levels of the ground state) and comparing the results. A Hilger 
medium quartz spectrograph was used for these measurements. 


§4. EXPERIMENTAL RESULTS 


4.1. The Ammonia—Oxygen Flame 
Ammonia and oxygen were reacted in the flat diffusion flame. ‘The appearance 
of this flame was described in our previous paper and two facts established by 
the earlier work were assumed in the present investigation, viz. (1) the emission 
in the flame is of thermal origin as proved by the fact that line reversal measurements 
with OH radicals and sodium atoms gave the same temperature, and (ii) the 


8 H. G. Wolfhard and W. G. Parker 


ammonia does not burn directly with oxygen but first undergoes thermal 
decomposition via NH, and NH. Quantitative measurements of the OH and O, 
molecules were made in the way described above. The results are shown 
graphically in fig. 3, which represents a cross section through the reaction zone 
1 cm above the top of the burner. The ammonia is on the left side and the oxygen 
on the right. ‘The temperature as measured by sodium line reversal method is 
also shown. 

Fig. 4(a) (Plate) shows a spectrogram of the visible and near ultra-violet 
region of an ammonia-oxygen flat diffusion flame. It should be read with the 
end view of the flame in mind (figs. 1 and 2). The abscissa is the plot across the 
reaction zone and the ordinate is the spectral wavelength. Beginning from the 
fuel side (on the left) the ammonia-« bands which are attributed to NH, radicals 
can be seen first, closely followed by the NH bands. After a clear gap in emission 
the main reaction zone follows with a faint continuum in the visible region caused 
by electron-ion recombination processes, and strong OH bands. Some bands 
of molecular oxygen can be seen a Jittle further towards the oxygen side. 

The ammonia absorption spectrum could not be measured quantitatively 
because it depends strongly on temperature. It was possible however to decide 
where all the ammonia was used up. This is denoted by the point A (fig. 3). 
As the ammonia approached the hot part of the flame it was observed in regions 
of fairly high temperature (>1000°c). Jt began to disappear rapidly at about 
2100°c at which point NH, and NH radicals appeared in emission. ‘This 
dehydrogenation was very unlikely to be a thermal dissociation in chemical 
equilibrium for several reasons. The equilibrium partial pressures of ammonia 
are negligible at temperatures of 1000°c and above even in an atmosphere of 
nitrogen and hydrogen. Moreover it is believed that the amount of NH (and 
probably of NH,) which accompanied the disappearance of the ammonia in the 
temperature region of 2100°c was also greater than the chemical equilibrium 
value. Although the NH bands are of thermal origin they are only strongly 
emitted in the position indicated in fig. 3. Some NH molecules may be present 
in equilibrium amounts all over the reaction zone and their existence is favoured 
in the reducing gases on the left side of the flame, but the amounts are exceedingly 
small as shown by Ward and Hussey (1949). The estimated partial pressure is 
of the order of 10~’ which is too low a concentration for detection by absorption. 
The non-equilibrium decomposition of the ammonia gas is thus established and 
is an important difference in the course of the reaction in this part of the flame 
compared with the main reaction zone where the hydrogen and oxygen react 
and chemical equilibrium is maintained. 

How the decomposition of the ammonia occurs is a matter for speculation. 
It is known that ammonia is very stable in the absence of surfaces and other 
catalysts and does not decompose with appreciable velocity below about 2000°c 
(Frank and Reichard 1936). In this flame dehydrogenation may be brought 
about by radicals diffusing from the main reaction zone which is some distance 
to the right as indicated by the peak of the temperature curve, fig. 3. It has 
been pointed out previously that the dehydrogenation does not involve simul- 
taneous oxidation to water since the oxidation process is strongly exothermic and 
the temperature distribution would differ from that found experimentally. ‘The 
actual heat liberation occurs further to the right (fig. 3). The decomposition of 
the ammonia therefore goes to N, and H, via NH, and NH. 
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The hydroxyl radical concentration was measured quantitatively from the 
absorption spectra in the way described in §3. This is expressed as the percentage 
partial pressure of OH (fig. 3). The maximum concentration is on the right 
of the position of maximum temperature, and this is to be expected since water 
(the main product of combustion) and oxygen in equilibrium give a higher 
partial pressure of OH than the combination of water and hydrogen which exists 
on the left of this temperature line. The strength of the OH emission bands 
depends on radical concentration and temperature and the peak is therefore 
shifted towards the position of maximum temperature. 

The oxygen concentration was measured quantitatively with two lines of 
different wavelengths. Below about 2% partial pressure of oxygen the absorption 
lines were too weak for measurement. Above about 40% of oxygen a difficulty 
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Fig. 3. Ammonia-—oxygen flat diffusion flame. 


arose from the nature of the flame. Absorption which is temperature-dependent 
occurred in a region of rapidly falling temperature on the outside of the flame. 
Reliable estimates of the oxygen concentration could not be made under these 
conditions. Between the limits of 2°% and 40%, the results with the two lines 
agree quite well, and at the position of maximum temperature the partial pressure 
was found to be about 4%. 

These results provide measured values of the temperature, and the concen- 
trations of OH radicals and oxygen molecules in certain parts of the main reaction 
zone. Although we do not know a priori the n,:n, ratio across the zone we 
can calculate the partial pressure of every component at any point in this zone 
provided our conception of the diffusion flame is correct, i.e. provided chemical 
equilibrium is effectively maintained in this part of the flame. ‘To calculate the 
concentrations of each component we need the temperature and the experimental 
value for the concentration of either the OH or oxygen molecules. Since the 
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results contain measured values of both these types of molecules in various parts 
it is possible to cross-check this method by assuming one component and calcu- 
lating the other. ‘This was done, and the partial pressures of OH radicals calculated 
from the measured values of oxygen are shown as full point circles in fig. 3. 
Agreement is good and may be fortuitous to some extent considering the possible 
errors involved in the measurement of the temperature as well as in the 
concentrations of the molecules, but it does suggest that the assumptions about 
the nature of the high temperature zone are correct as regards the oxygen diffusion 
flame. 

The partial pressure of the oxygen could not be determined experimentally 
below about 2% for the reasons given above. OH concentration was still high 
in that part of the flame however and the oxygen curve was therefore extended 
(dotted portion of fig. 3) by calculation from the measured temperature and 
OH concentration. The partial pressure of molecular oxygen at the point where 
the ammonia decomposes and the NH absorption bands are strongest is thus 
found to be about 0:01%. 

We conclude from these results that the ammonia—oxygen diffusion flame 
consists of two parts (i) a region of thermal breakdown of ammonia to nitrogen 
and hydrogen and (ii) a high temperature region where the hydrogen burns in 
oxygen and nitrogen is a diluent. Chemical equilibrium is virtually maintained 
in this region but a small chemical potential exists in (i) where the temperature is 
above the equilibrium decomposition temperature of the ammonia. ‘This occurs 
outside the main reaction zone of the flame and the chemical potential involved 
cannot be compared with that in the inner cone of a pre-mixed gas flame. 

One other observation in the flame is worth mentioning. ‘The oxygen bands 
most readily obtained in emission were (0, 13), (0, 14) and (0, 15) and these bands 
did not absorb appreciably because the high levels in the ground state to which 
the transition occurs were not appreciably excited. It was therefore possible to 
calculate the relative strength of the emission bands (assuming that they were 
of thermal origin) since the temperature and the oxygen concentration had been 
measured and the amount of oxygen in the vibrational ground level of the n° = 
state was therefore known. These calculated values are inserted in the graph 
(fig. 3) and have been given an arbitrary maximum. The position of the maximum 
is fixed by the relative effects of oxygen concentration and temperature. The 
graph also shows the experimental peak for the strength of the O, bands denoted 
by B. Exact measurements of the O, emission were rather difficult because the 
bands were onastrong background. The positions of the peaks of the experimental 
and calculated values are sufficiently close together, however, to suggest that the 
oxygen is really thermally excited; the implication of this is discussed later. 


4.2. The Ethylene-Oxygen Flame 


‘The ethylene-oxygen flat diffusion flame was similar to the methane-oxygen 
flame examined in our initial work and had a number of advantages for the 
present investigation. ‘I'he reaction zone was thicker and very stable and carbon 
formation was more pronounced and could be readily observed in absorption. 
An end view of the flame is shown in fig. 2. 

Fig. 4(b) shows the emission spectra of an ethylene-oxygen flat diffusion 
flame. ‘This too should be read with the end-view of the flame in mind. 
Beginning on the left, 1.e. on the fuel side of the flame, one can see first the carbon 
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continuum and closely attached to that are the C, and CH bands. NH bands 
which are due to impurities are strongest in the same position as the CH bands 
but extend well into the main reaction zone. Between the carbon continuum 
and the main reaction zone there is a well marked dark space. ‘The main reaction 
zone gives principally a broad continuum which has its origin presumably in 
an electron-ion recombination process. The OH bands are most prominent 
in this part but it is remarkable that some OH bands can also be seen right over 
on the fuel side. ‘he oxygen bands may be seen on the oxygen side of the flame. 

The temperatures in this flame were higher than those of the ammonia flame 
and determination of the temperature curve required the use of a carbon arc. 
‘The carbon arc was first calibrated against a standard tungsten lamp using a 
stepped sector. Flame temperatures determined by OH-line reversal and Na-line 
reversal were again identical indicating thermal equilibrium. OH and O, 
concentrations were measured as in the case of the ammonia flame and the results 
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Fig 5. Ethylene-oxygen flat diffusion flame. 


are plotted in fig. 5. One of the most striking things about these results is the 
large distance between the point at which the O, concentration becomes very 
small and the position of the luminous carbon zone. This means that little or 
no oxygen reaches a part of the flame where the hydrocarbon undergoes extensive 
chemical change. This can also be seen in fig. 4(c). ‘The oxygen absorption 
lines fail to reach the region marked on the spectrogram which represents the 
beginning of the absorption by the carbon zone. 

Again it is believed that the hot gases in the main reaction zone are in chemical 
equilibrium. The OH partial pressures were not calculated from the measured 
oxygen concentrations as a method of cross-checking the experimental results 
in this system for two reasons: (a) The calculation requires a very exact knowledge 
of the temperatures in the zone and the experimental temperature curve is not 
so accurate as in the ammonia—oxygen flame; this is because the carbon arc had 
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to be used as the background source for the line reversal method. (b) An unknown 
proportion of the hydrocarbon is precipitated as carbon in the carbon zone. 
This carbon is effectively lost and the exact 7,: my ratio is not known on which to 
base the calculated chemical equilibrium. 

The region in which carbon formed in the flame was examined in more detail. 
The absorption of the carbon zone was very strong. ‘I'he absorption occurred 
somewhat on the left of the emission because of the marked temperature gradient 
in that part of the flame. ‘The width of the absorption, however, was fairly big 
and it is doubtful whether the whole of this absorption was due to carbon particles. 
On the left side of the absorption the temperature is still high and about 2000°c 
according to an extrapolation of the temperature curve in fig. 5, yet no continuous 
radiation was observed. It may also be noted that in the position of this absorption. 
there was a reappearance of OH emission bands. The OH bands, which were 
strongest in the reaction zone, had faded away before reaching the carbon zone 
(fig. 4(d)). Their reappearance on the fuel side is a curious effect which only 
occurred when there was strong carbon formation. It is thus unlikely to be a 
true radiation but scattered light which has its origin in the main reaction zone. 
The nature of the particles responsible for scattering is a matter of conjecture. 
Their size must be less than a wavelength. They might be carbon particles in 
which case one would expect them to be red hot and radiate at the prevailing 
temperature or they might be pre-carbon particles with small emissivity. The 
latter would account for the fact that absorption occurs 1 mm on the fuel side of 
the radiant carbon zone. Such particles will scatter the radiation of the background 
source and produce an apparent absorption which is not a true absorption obeying 
Kirchhoft’s law. 

The width of the absorbing carbon zone gradually became much greater 
below 25004, an effect which also occurred in a flame of methane saturated with 
benzene vapour and may be seen in the spectrogram, fig. 4(e). This certainly 
represents an additional absorption effect which is not directly connected with 
the carbon zone. ‘This additional absorption has a wedge-shaped appearance 
and seems to be related to the absorption previously observed by the authors 
in pyrolysis experiments on pure hydrocarbons (Parker and Wolfhard 1950). 
It is therefore proposed to call this additional absorption ‘ pyrolysis absorption’ 
in contrast to the absorption in the carbon zone. Straight chain and aromatic 
hydrocarbon gases showed such an absorption on heating which spread towards 
longer wavelengths at higher temperatures. This phenomenon was attributed 
to intermediate compounds which lead to carbon formation. 

It is well known that hydrocarbons alone on heating to 1 000° c or over produce 
carbon, and the conditions in the diffusion flame in the region considered here 
are mainly thermal since no oxygen reaches it. Water vapour, carbon dioxide 
and carbon monoxide are present however and diffuse from the main reaction 
zone. ‘I'hese gases may react with the carbon particles but their influence as 
far as pyrolysis and carbon formation are concerned is probably small and largely 
a dilution effect. We therefore suggest that the carbon formation in an oxygen 
diffusion flame occurs in the same way as it does in the pyrolysis of a hydrocarbon. 

Analysis of spectra was limited on the fuel side of the flame because ethylene 
does not absorb in the visible and ultra-violet regions. It was not possible 
therefore to follow its thermal changes in the same way as with ammonia. Acetylene 
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and more especially benzene respond in these regions and, since their reactions 
in the diffusion flame are probably similar, a few experiments were made with these 
compounds. Both produce large and embarrassing quantities of carbon in the flat 
diffusion flame so that acetylene was used diluted with hydrogen, and benzene was 
introduced into a methane flame by passing methane gas through a benzene 
saturator. 

The absorption bands of acetylene were observed on the fuel side of the 
flame and reached the region of continuous absorption described above as 
pyrolysis absorption. Benzene, on the other hand, which absorbs strongly at 
wavelengths of the order of 26004 is not obscured by the pyrolysis absorption 
as the latter occurs at shorter wavelengths. There was a gap between the break- 
down of the benzene as indicated by the disappearance of its absorption, and the 
onset of the continuous absorption of the carbon zone. It seems clear from these 
results that the hydrocarbon has ceased to exist shortly after the commencement 
of the pyrolysis continuum. The distance which the undecomposed benzene 
diffuses can be seen relative to the position of the carbon zone in the spectrogram 
of the flame (fig. 4(e)). The reactions of the hydrocarbon are thermal (since 
measurements already described indicate that free oxygen and OH radicals are 
not present in this part of the flame) and lead to hydrogen and solid carbon in 
the radiant zone at 2300°c. Not all the hydrocarbon is degraded into carbon 
however, otherwise methane would give about half the amount of carbon that 
comes from ethylene or acetylene; the observed ratio is more like 1:10. 

It is possible that the stationary concentration of solid carbon is controlled by 
the water vapour present and that the decomposition of the hydrocarbon into 
solid carbon and hydrogen is complete in the carbon zone at 2000°c to 2 300° c 
but the carbon reacts with the water vapour to form hydrogen and carbon monoxide 
which diffuse away and are burned in the main reaction zone further to the right. 
This process might be relatively more important in the case of methane than, 
say, acetylene because the water vapour concentration is obviously higher. 

As soon as carbon particles have formed, their size prevents their transport 
horizontally by diffusion. The radiant carbon zone beginning near the burner 
top represents therefore a stream line in the flame and all other spectroscopic 
phenomena may conveniently be considered in relation to the position of the carbon 
zone. 

A little to the right of the radiant carbon zone (i.e. on the oxygen side) rather 
weak C, bands occurred. These bands appear strongly in flames of pre-mixed 
hydrocarbon air or oxygen gases but their origin is not understood (Gaydon 
and Wolfhard 1950). It has been suggested that polymerization of the C, 
molecules is the source of carbon in the flame. Apart from energy considerations 
and other objections to this suggestion (Parker and Wolfhard 1950) the situation 
of these bands on the right of the carbon zone proves this suggestion is wrong. 
Both carbon and C, molecules need the same amount of energy for their excitation 
at the same wavelengths and we should expect the C, bands to appear inside the 
carbon zone or even a little to the left if C, molecules are formed before solid 
carbon. The different intensities of emission of C, in pre-mixed and diffusion 
flames is also worth noting. In pre-mixed flames the C, emission gets its strength 
from non-equilibrium conditions in the flame; in diffusion flames the bands 


are much weaker. 
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It seems to the authors that a more likely explanation at least of part of the 
C, bands is that C, molecules in the ground state are produced by vaporization 
of the carbon particles and that these molecules are thermally excited afterwards. 
It is known that C, radiation may be observed at corresponding temperatures 
in carbon furnaces (Brewer et al. 1948). The number of C, molecules does 
not therefore depend on the number of carbon particles but only on the presence 
of carbon and the temperature. C, will thus radiate most strongly at the right- 
hand side of the carbon cloud and the small shift of the C, bands to the right of 
the radiant carbon zone is readily explained on account of the temperature 
gradient. Two other observations are also accounted for. Firstly the C, bands 
were weaker in the methane flame than in the hotter ethylene or acetylene flames 
and secondly the introduction of small quantities of benzene into the methane 
flame did not affect the intensity of the C, bands. In the latter case the amount | 
of carbon increased considerably but the temperature remained about the same. 

It is still doubtful, however, whether all the C, radiation can be explained 
in this way in spite of its relatively feeble strength. Observations at different 
heights in the flame showed that the C, and CH radiations were stronger at the 
base of the flame than at the level normally investigated 1 cm above the burner. 
This was no surprise since some pre-mixing may be assumed in the ‘dead space’ 
at the burner mouth, and the characteristics of a pre-mixed flame in which C, 
and CH are fairly strong was to be expected. The reduction in strength of the 
C, and CH bands was rather gradual however in the vertical direction above the 
‘dead space’. It is therefore possible that more than one process of C, formation 
and excitation is involved. A distinction between the physical process of 
evaporation and thermal excitation of C, and the chemical process of C, formation 
similar to that in the pre-mixed flame might be attempted by examining the 
rotational temperature. It has been shown by Gaydon and Wolfhard (1950) 
that the rotational temperature of the C, is abnormally high in the pre-mixed 
gas flame. ‘The application of this test would be difficult however since the 
C, radiation is very weak and close to the strongly radiant carbon zone. 

‘The reason for the intensity of the CH radiation is also obscure. The absence 
of CH bands from the spectra of the interconal gases of pre-mixed gas flames 
indicates that the probable chemical equilibrium concentrations of this radical 
(Ward and Hussey 1949) are too small to produce any radiation. CH bands 
are observed in the inner cones of pre-mixed flames but they are not of thermal 
origin (Gaydon and Wolfhard 1948). ‘The rotational temperature estimated 
from them is not abnormal and there is consequently no easy method of dis- 
tinguishing between non-equilibrium and thermal radiation. In the diffusion 
flame the CH radiation was too weak for temperature measurements by the 
line reversal method and it was therefore impossible to obtain a value for 
comparison with the known temperature in that part of the flame. We suggest 
that the CH is formed chemically by a process such as C,+OH-~CH+CO. 
This reaction is probably very exothermic and may produce CH radicals which 
are chemiluminescent. On the other hand the CH radicals formed by this 
process may be thermally excited subsequently. More experimental evidence 
is required to support or refute the suggestion. Im the present work we noticed 
that the strength of the CH radiation varied at different heights in the flame. 
It behaved in a similar manner to the C, radiation and was stronger at the base 
than it was higher up the flame. 
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4.3. Additional Observations from the Spectra of other Flat Diffusion Flames 


(a) Hydrocarbon-nitric oxide flames. A hydrocarbon flat diffusion flame could 
be burnt with nitric oxide instead of oxygen. It seems that combustion occurred 
because the nitric oxide had time to dissociate extensively into nitrogen and 
oxygen, and combustion in the main reaction zone was essentially similar to that 
in the oxygen flame. ‘This could be proved by the presence of lines of the oxygen 
bands in the main reaction zone. The end view of the flame was also similar to 
the oxygen flame and apparently divided into the main reaction zone (blue 
continuum) and a carbon zone (white light). Decomposition of nitric oxide 
requires a high activation energy and decomposition was not complete even in 
the main reaction zones. - Nitric oxide in quantities greater than equilibrium 
persisted and diffused into the carbon zone. By using a methane-nitric oxide 
flame which had only a moderately intense carbon zone it was possible to measure 
the absorption bands of nitric oxide in the luminous carbon zone. Other bands 
were also observed and compared with similar bands in the oxygen flame. The 
CH emission bands had not changed but NH bands were of course stronger in 
the nitric oxide flame. In the oxygen flame these bands are due to a small amount 
of nitrogen from impurities. ‘The strength of the carbon emission was perhaps 
a little greater in the NO flame and this point was carefully noted because nitric 
oxide is known to act as a chain-breaking molecule. Introduction of nitric oxide 
as diluent into the fuel of a methane—oxygen flame gave essentially the same 
result. If carbon formation is a chain reaction the nitric oxide might be expected 
to cause a reduction in carbon formation when oxygen is replaced by nitric oxide 
or when nitric oxide is introduced into the fuel. This effect was not observed. 


(b) Carbon monoxide—oxygen diffusion flame. ‘The flat diffusion flame of carbon 
monoxide gave principally radiation in the main reaction zone and this consisted 
of a strong continuum extending far into the ultra-violet region. At about 
40004 there were some bands superimposed on the continuum and these were 
identified as the so-called carbon monoxide flame spectrum. ‘The bands are 
probably emitted by CO,. Besides this system there were some OH bands due 
to moisture and these could be seen even in absorption. ‘The Schumann—Runge 
oxygen bands were also present. ‘The relative positions of these bands is of 
interest because Hornbeck has recently suggested that part of the reported CO, 
band structure is attributable to the O, bands (Hornbeck and Hopfield 1949). 
Fig. 6, which records the results obtained from the spectrogram of the carbon 
monoxide-oxygen flame, indicates that this suggestion is wrong because the two 
systems may be seen separated in space. ‘The oxygen bands occur much more 
on the oxygen side than the CO, bands and their separation in this way 1s proof 
of their independence. 


(c) The hydrogen-oxygen diffusion flame. ‘This flame was considered in the 
previous paper and it was mentioned that the oxygen bands appeared strongly 
in emission and could be observed from 4000 A down to 2400 and even lower. 
It has since been suggested that excitation of oxygen molecules in flames is not 
of thermal origin because the Schumann—Runge bands require a very high 
excitation energy. Thus in the carbon monoxide—oxygen pre-mixed gas flame 
Laidler (1949) proposed the following mechanism: CO,* + O, > O,*+CO, in 
which an excited CO, molecule hands on its energy during collision to an oxygen 
molecule. The hydrogen flame provided an opportunity for deciding whether 
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the O, radiation was thermal er due to chemiluminescence, because the oxygen 
bands stood out clearly on the spectrogram and were not partially obscured by 
a continuum as in the carbon monoxide flame. The method employed was to 
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Fig. 6. Relative position of OH, O, and CO flame bands in the carbon monoxide-oxygen ) 
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measure the excitation temperature by the line reversal technique and compare 
the result with the theoretical flame temperature. A hydrogen discharge lamp 
was used as the background light source and this was previously calibrated with 
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a carbon arc. ‘The brightness of the carbon arc was assumed to be 3800°xK 
independent of wavelength. ‘The line reversal measurement was not easy because 
it required a long exposure time. The accuracy was probably not high but a 
reversal temperature of 2 800°c + 60°c was estimated and this agreed well with 
the temperature of the main reaction zone. The oxygen molecules like the 
OH molecules were therefore thermally excited in the flame, a result which was 
also suggested by the position of the O, emission bands in the ammonia, and 
ethylene—oxygen diffusion flames. The oxygen bands in the carbon monoxide- 
oxygen flame were weaker than those appearing in the hydrogen—oxygen flame. 
Hence there is no reason to suppose that the oxygen molecules are not in 
thermal equilibrium in the CO flame as well. 


§5. A COMPARISON OF PRE-MIXED AND DIFFUSION FLAMES 

The experimental results which have been obtained with oxygen diffusion 
flames illustrate some of the fundamental differences between these flames and 
those of the Bunsen type. 

Firstly all pre-mixed gas flames start with a high chemical potential which 
breaks down in the reaction zone. Beyond that zone, i.e. above the inner cone 
of the Bunsen flame, the hot gases are essentially in chemical equilibrium. The 
diffusion flame on the other hand never develops large chemical potentials in the 
main reaction zone. Diffusion endeavours to change the gas composition but 
chemical reaction occurs at once and restores the equilibrium. ‘This is probably 
only true in the main reaction zone and some chemical potential may develop 
during the decomposition of the fuel and oxidizer in the preheating zones on either 
side of the main reaction zone. Non-equilibrium conditions were observed in 
the decomposition of ammonia as it approached the main reaction zone and 
probably also in the decomposition of nitric oxide. 

Secondly, there was no evidence of chemiluminescence in the main reaction 
zone of the diffusion flame, and no evidence of intermediate radicals such as HO, 
or HCO was found in either emission or absorption spectra. ‘The absence of 
HCO bands in the main reaction zone of the hydrocarbon diffusion flame is 
understandable, since the principal reactions taking place there are between 
hydrogen, carbon monoxide and oxygen. 

The reaction in the main reaction zone is a relatively slow process. If 
molecules such as OH, H,O, CO, etc. are formed in an excited state by the 
reaction they are unlikely to be detected in that state because of the much stronger 
thermal radiation of the majority of the molecules which are already present in 
chemical equilibrium. The chemiluminescence of molecules such as CH, 
HO, etc., if formed as intermediates, should be discernible, however, because 
these molecules are never present in equilibrium. If the HO, radical is a 
comparatively long-lived intermediate of the combustion of hydrogen in oxygen 
it should be detected by thermal radiation or absorption in a flame of this kind 
where the long optical depth makes the test very sensitive. No absorption or 
emission which could be attributed to HO, was in fact found. It was characteristic 
df the diffusion flame that the C, and CH radiations were not observed in the main 
reaction zone, but occurred outside and in close proximity to the carbon zone. 

Chemiluminescence is commonly observed in pre-mixed gas flames though it 
is not the only reason for enhanced radiation. Pre-mixed gas flames of hydro- 
carbons have a high excitation temperature which shows up particularly in the 
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ultra-violet where the thermal radiation diminishes. This effect may be demon- 
strated by introducing iron atoms through the medium of iron pentacarbonyl 
(Gaydon and Wolfhard 1951). The intensity ratio between the iron lines in 
the far ultra-violet and those in the visible region is abnormal and is comparable 
with the corresponding ratio obtained in the spectrum of an iron arc. Some 
lines are more than a thousand times as strong as they should be. In diffusion 
flames at atmospheric pressure the results were very different. The iron behaved 
as though it was in thermal equilibrium and the iron lines in the far ultra-violet 
were very faint. This is an important difference between the two types of flames. 
It probably arises because in pre-mixed gas flames intermediates such as peroxides 
occur which cannot be formed in diffusion flames because the oxygen does not 
come in direct contact with the hydrocarbon molecule. 

This fact may also explain the thermal behaviour of the OH radiation in 
diffusion flames, which again differs from that of pre-mixed gas flames. In the 
latter the OH radiation is chemiluminescent and possesses a high rotational 
temperature. It may originate in a reaction such as CH+O,>CO+OH*. 
This reaction is very improbable in a diffusion flame since the CH radiation is 
relatively weak and occurs at a position in the flame where the oxygen concen- 
tration is already very low. The OH radiation is thus confined to a region where 
the OH concentration and the temperature are high. ‘These conditions exist 
more on the oxygen side of the flame. 

Diffusion flames of hydrocarbons and oxygen showed certain changes at very 
low pressure. In contrast to the pre-mixed gas flame the thickness of the reaction 
zone did not alter with the pressure but the main reaction zone and the hydrocarbon 
decomposition zone appeared to merge slowly into one another. Carbon formation 
became weaker and ceased altogether at very low pressure. C,, CH and OH 
radiations grew stronger. Iron pentacarbonyl! was introduced and the iron lines 
again showed a high intensity in the far ultra-violet region indicating a high 
excitation temperature at low pressures. ‘This change in behaviour of the 
diffusion flame at low pressure fits into our conception of the flame (§2 ef seq.). 
At low pressure the reaction velocity decreases while the diffusion rate increases. 
The gases in the flame will begin to mix and allow a certain chemical potential 
to develop before the reaction can complete. The characteristics of the flame wiil 
change and some of the characteristics of the pre-mixed gas flame will appear. 
‘The increased strength of the C, and CH radiations, the high rotational temperature 
of the OH molecules, and the high excitation temperature of the iron atoms amply 
support this theory. 

Lastly, a clear difference between diffusion and pre-mixed gas flames is shown 
by the influence of the oxides of sulphur on carbon formation. In pre-mixed 
gas flames sulphur trioxide promotes carbon formation whereas it has the reverse 
effect in diffusion flames (Wolfhard and Parker 1950). 


§6. CONCLUSION 
Visual observation of the flat flame immediately reveals two distinct zones. 
In the ammonia-oxygen flame there is broadly speaking a yellow zone on the fuel 
side and a blue zone on the oxygen side with a thin dark space between them. 
In the hydrocarbon-oxygen flame there is also a yellow zone on the fuel side 
which is easily recognized as the emission of incandescent carbon particles, and 
again there is a blue zone on the oxygen side with a dark space between. The 


Structure of Diffusion Flames 19 


similarity between the blue zones in the two flames is very marked and these 
zones resemble the single zone observed in the hydrogen-oxygen diffusion flame. 
Optical measurements revealed a great deal of more precise detail and it has 
now become necessary to label the different regions of the flame in order to avoid 
_ confusion. The blue zone is better identified as the main reaction zone because 
it is the hottest part of the flame. The emission from this zone is thermal and 
the constituent gases are essentially in chemical equilibrium. On either side of 
_ the main reaction zone a steep temperature gradient exists and these outer regions 
are broadly termed pre-heating zones. Fuels such as hydrocarbons and ammonia 
and oxidizers such as nitric oxide undergo considerable chemical change in these 
pre-heating zones. On the fuel side the change is principally a thermal 
decomposition because the main reaction zone presents a barrier to oxygen 
"penetration. 
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ABSTRACT. Processes of electronic excitation in flames are briefly discussed, and it is 
shown that if there is lack of equilibrium the reversal temperature may be closer to the 
effective vibrational temperature of the surrounding molecules than to the mean translational 
temperature. The lack of radiative equilibrium in all flames will cause measured tempera- 
tures to be too low; for a flame with air at 2 000° k using Na the error will be between 
3-3° and 1-6° at 1 atm. and more at lower pressure. For hot flames, above 2 700° k, the 
OH band may be used for reversal measurements; it appears to give satisfactory results 
and has some advantages over the use of sodium. 
2-2 
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§1. INTRODUCTION 

HE spectrum-line reversal method of measuring flame temperature is 

experimentally quite simple and in fairly common use. The method 

depends on Kirchhoff’s law that in equilibrium the emissivity and 
absorption coefficient at any wavelength must be equal. A method depending 
on this law was first used for luminous flames containing carbon particles by 
Kurlbaum (1902), and in 1903 Féry first used the resonance lines of sodium and 
other metals. The spectrum of a controllable continuous source, such as a 
strip-filament tungsten lamp, is viewed through the flame to which a little sodium 
or other metal has been added, and the brightness of the background source is 
adjusted so that the spectrum line of the added metal is seen neither in emission 
nor absorption by the flame against the background; the background source 
is then assumed to have the same brightness temperature as the flame, and hence 
the flame temperature can be determined by calibrating the background source 
against a black body. The experimental method has been fully described in 
a number of papers (Griffiths and Awbery 1929, Henning and Tingwaldt 1928, 
Loomis and Perrott 1928). 

The method depends on the existence of equilibrium in the flame gases, and 
for flames with air numerous comparisons have been made between temperatures 
determined by the spectrum-line reversal method and those determined by 
other methods or calculated from thermochemical data. In the interconal gases 
of a premixed flame, well above the reaction zone, agreement is generally 
satisfactory and there is little reason to distrust the method. For the actual 
reaction zone of a flame, equilibrium is not to be expected, and it would not be 
surprising if the method gave results quite different from those expected for 
theoretical flame temperatures; indeed it is not easy to define temperature in 
the reaction zone of a flame. Our measurements (Gaydon and Wolfhard 1951) 
in the reaction zones of low-pressure flames have shown quantitatively the 
existence of abnormally high electronic excitation, and we have also shown more 
qualitatively that such effects also exist in the reaction zone at atmospheric 
pressure. 

One of the most interesting regions of a flame is that just above the reaction 
zone or inner cone. Here disagreement between theoretically expected 
temperatures and those measured by spectrum-line reversal have been reported. 
Thus Lewis and von Elbe (1943) found that for flames with air the measured 
temperature was appreciably below the theoretical temperature in the region 
just above the visible reaction zone, and did not attain the theoretical value 
until some 6 to 10 mm above the inner cone. On the other hand, for flames 
of lean mixtures with oxygen, the temperature in this region appears too high. 
The reason for these discrepancies is not known, although Lewis and von Elbe 
suggest that it may be due to lag in vibrational excitation of the molecules. 

With sodium, absorption in the cooler outer layers of the flame is always 
likely to give too low a temperature unless precautions are taken. Strong, 
Bundy and Larson (1949) have made an experimental study of this effect using 
an interferometer to examine the contours of the spectrum lines. 

Experimental work on the sodium-line reversal method for flames with air 
indicates that provided we are well away from the reaction zone, the method 
gives reliable results; for very hot flames with oxygen not much work has been 
done to test it. The physical processes causing electronic excitation in flames 
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do not seem to have been very fully considered. It is our purpose here to discuss 
these, and the effect that departure from true equilibrium might have on measured 
temperatures by the reversal method. We shall also discuss the magnitude of 
errors due to the flame gases not being in radiative equilibrium. For very hot 
flames there are some practical difficulties in using sodium, and we shall indicate 
the advantages we have found in using lines of the OH band for the reversal 
method. 


§2. MECHANISM OF ELECTRONIC EXCITATION 

In electric discharges and arcs the excitation to higher electronic states is: 
caused by impact of free electrons accelerated in the electric field. In flame 
gases, however, the concentration of electrons is much lower and there is no: 
field. ‘The number of collisions with ordinary molecules is very much greater 
than with electrons (by a factor of at least 10°), and since we shall see that 
excitation by molecules may be fairly efficient, that due to free electrons must 
be negligible in comparison. 

Conversion of energy of electronic excitation into other forms can be studied 
by the quenching of fluorescence, and then, in equilibrium, we may apply the 
principle of microreversibility to learn about the reverse processes causing 
electronic excitation. 

We should not expect that large amounts of translational energy could readily 
be converted into that of electronic excitation, and this is confirmed experimentally 
by the extreme weakness of quenching by monatomic gases. On the other hand, 
molecules, either diatomic or polyatomic, are frequently quite efficient in quenching 
fluorescence. The efficiency varies a great deal with the nature of both the 
quenching molecule and the excited species. The efficiency may be expressed by 
the value of the collision cross section o,,,”._ "The number of collisions a molecule 
of type A suffers per second with molecules of type B is 


F 8nRT(M, +My) |??? 
FaBalte asa | 
where m,, is the number of molecules of B per cm’, R is the gas constant, T the 
absolute temperature and M, and M,, the molecular weights of A and B. 
Laidler and Shuler (1951) have collected some data (from Norrish and Smith 
and from Terenin and Prileshajewa) for quenching cross sections in various 


gases. For excited sodium Na(2P) the following values are given, in (Angstréms)? : 


Inert gases 0 CH; 0-11 
H. 74 C,H, 0-17 
Nz 14-5 C,H, 0-2 
CO 28-0 C,H, 44 

\e 40-0 C,H; 52 
Br, 100-0 


Cross sections for elastic collisions from viscosity data are of the order 10A?, 
varying of course with the type of molecule. For an excited alkali atom they 
would be relatively large. Thus we see that the quenching efficiency may be of 
the order unity, or may be quite small. The general mechanism of energy 
transfer on collision has been reviewed by Laidler and Shuler (1951) and by 
Massey (1949). The results are best interpreted by considering the potential 
energy surfaces of the temporary collision complex. If there is an intersection 
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of surfaces between which a radiationless transition is probable, then there is 
a chance that during the collision such a transition will occur, and as the complex 
breaks up it may then leave the products in different electronic states. ‘This 
may lead in some cases to a dissociation of the colliding molecule, if the energy 
is sufficient, but more often the energy set free by the change from a high to a 
lower electronic state will be released mainly in the form of internal vibrational 
energy. We can see that the efficiency of such processes will depend on the form 
and position of intersection of the potential surfaces and may therefore vary 
widely with the nature of the colliding particles. 

We may summarize this section by saying that electronic excitation in flames 
is unlikely to be due to collision with free electrons, or to conversion of trans- 
lational energy into energy of electronic excitation. It may result from chemical 
processes, such as recombination of free atoms or radicals, and this type of process 
is likely to be especially important in the reaction zone, but the main process 
responsible for maintenance of equilibrium between excited and unexcited atoms 
in the interconal gases will be conversion between vibrational and electronic 
energy. 

Thus if there are small departures from equilibrium, due perhaps to lag in 
vibrational excitation, then we should rather expect that the temperature given 
by the spectrum-line reversal method may be nearer to the mean effective 
vibrational temperature of the surrounding molecules than to their mean 
translational temperature. A measured temperature above the theoretical 
maximum value may therefore indicate that the molecules have some excess of 
vibrational energy and a deficit of translational energy. 


§3. EFFECT. OF DEPARTURE FROM RADIATIVE EQUILIBRIUM 

Complete equilibrium is only attained in a constant temperature enclosure 
such as a black-body furnace. ‘The majority of laboratory flames are optically 
thin, and depopulation of excited electronic states by emission of radiation is 
not compensated by absorption of radiation. The spectrum-line reversal method 
essentially determines the ratio of the populations of the excited and ground 
electronic states, and any reduction in the population of the excited state by 
radiation will cause the measured temperature to be lower than the mean 
temperature of the gas molecules. 

Now we may assume that the main process responsible for maintaining 
equilibrium is the collision with gas molecules. In equilibrium the number of 
atoms excited and de-activated by collision must be equal. If we know the 
quenching cross section we can calculate the average time before de-activation 
by collision, and we can then compare this with the average time before de-activa- 
tion by radiation if we know the radiative lifetime. ‘To obtain an estimate 
of the magnitude of such effects, we shall consider sodium atoms in a flame with 
air at Z2000°, and assume that the majority of colliding molecules are nitrogen. 

Taking the value of 14-5 x 10-16 cm® for the quenching cross section and using 
the formula for number of collisions, we find that the excited Na atom makes 
3-1 x 10° effective quenching collisions per second. For sodium the radiative 
lifetime (Stephenson 1951) is 1-6 x 10-8 sec. Thus the amount of depopulation 
by radiation, in the absence of any re-absorption, may amount to 2:0% of the 
depopulation by collision. Now from the Maxwell—Boltzmann law the pro- 
portion of atoms in the excited state should be exp(—E/kT). If we put 
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E=16964 cm™ for Na(??P), R=0-697 cm! and T=2 000, we obtain a value for 
the population in the excited state. If we then assume the 2:0°% reduction in 
population by loss of radiation and use the same expression to deduce the 
temperature, we find that it is 3-3° lower. A similar calculation for a flame at 
1 


jp atm. pressure, when the number of quenching collisions will be smaller, 
1 


gives an error through radiation loss of 27°, and at + atm. of 160°. 

These calculations have required a number of assumptions. First we assumed 
that the surrounding gas was nitrogen. It is possible that quenching by other 
molecules like H,O and CO, might be rather different, but since that for N, is 
around unity, the effect would be small. We have also tacitly assumed that the 
quenching cross section in a hot gas is the same as for quenching experiments 
in cooler gas; again, since the efficiency is near unity, any change is likely to be 
small. We have also assumed no self-absorption; now for the reversal method 
to work, absorption must be appreciable, at least in the centre of the line. Even 
for a completely opaque flame, the outer layers of gas, which are the layers which 
would then be effective, would receive radiation from only a hemisphere while 
radiating to a whole sphere, so about half the calculated effect would occur. 
Also the wings of the line will not be self-absorbed to any extent, and will tend 
to show the whole effect. ‘Thus the actual error will be between the calculated 
one (3-3° at 1 atm.) and a half of this value (1-6° at 1 atm.). 

It is interesting to have some quantitative knowledge of this error. For 
flames of different composition, or for different spectrum lines, the error would 
be rather different; it seems that it is not likely to be very serious at atmospheric 
pressure, but would tend to become more important at low pressure. 


§4. THE USE OF THE OH BAND FOR REVERSAL WORK 

While sodium is the most convenient metal to add to a flame with air, for very 
hot flames, such as those of near stoichiometric mixtures with oxygen, it is not 
easy to introduce a controlled amount of sodium. For these powerful flames 
there is always a risk of the flame striking back violently; this complicates 
the use of spraying devices, and the flames are too hot to introduce the sodium 
directly ona wire. Also the sodium-line method always has the disadvantage that 
the sodium in the cooler outer layers of the flame tends to contribute to the 
result. 

Now practically all flames contain OH radicals, and the OH band at 3064 A 
can be reversed. Generally the concentration of OH tends to be highest in the 
hottest part of the flame, so that the effect of the cooler outer layers is less 
important than with sodium. Also by using lines of high rotational energy one 
can exercise a further selection in favour of the high temperature region of the 
flame, since levels of high rotational energy will only be highly populated in the 
hottest part of the flame. 

For flames above 2 700° kit is not possible to use a tungsten lamp as background 
source, and a carbon arc is the most suitable source yet developed. ‘The positive 
crater has a brightness temperature at 6500A of about 3818°k (Henning and 
Tingwaldt 1928). Since the brightness of the arc itself cannot be varied, the 
effective brightness may be reduced with an adjustable rotating sector. Here 
again the OH reversal method is more sensitive than the Na method because 
at high temperatures a given change of brightness, say by 10%, corresponds in 
the yellow to a fairly large change in temperature, but in the ultra-violet around 
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3100 it corresponds to a much smaller change in temperature; thus if we can 
fix the reversal point to within an intensity of +10, the error in terms of 
temperature will be much smaller for OH than for Na. 

The OH reversal method has the disadvantage that it has to be done photo- 
graphically and with fairly large dispersion. We have found a quartz Littrow-type 
spectrograph, such as Hilger’s commercial instruments, quite suitable. A medium 
quartz spectrograph can be used, but is less sensitive. | With the Littrow-type 
instrument, exposure times for flames at around 3000°K are of the order 
10 seconds. For flames below 2700° exposure times would be very long, and 
the sector used in front of the arc would have to be set so finely that it would 
become difficult to read it accurately. The optical arrangements are of the usual 
type, except that more reliable results are obtained by using a quartz-fluorite 
achromat lens because the focus of the arc on the slit may be poor in the ultra-violet 
if it is adjusted visually and an achromat is not used. 

We have used the OH reversal method on a number of flames, and for lines 
of high rotational quantum number have always obtained temperatures for the 
interconal gases which are in agreement with theoretical temperatures within 
the limits of experimental error. For lines arising from levels of low rotational 
quantum number we obtain fair agreement, but there is a tendency for the values 
to be somewhat low, by up to 30°; this is probably due to absorption in the 
cooler outer layers of the flame. In and very near the reaction zone the reversal 
temperatures are found to be too high. 

The radiative life for excited OH is about 4 x 10~® sec, which is much longer 
than that for Na. We might thus expect the lack of radiative equilibrium in flame 
gases to have a smaller effect on OH; however, we have obtained some evidence 
(Gaydon and Wolfhard 1948) that the quenching of excited OH by flame gases 
has, in terms of a normal gas-kinetic collision cross section, an efficiency of only 
about ;, so that the longer effective collision life may cancel the good effect of 
the longer radiative life. ‘The OH is likely to be more affected than Na by the 
abnormally high electronic excitation in the reaction zone; we have shown that 
this becomes more important at shorter wavelengths. However, this effect does 
not appear to extend more than about twice the thickness of the reaction zone, 
and should therefore not affect the main body of interconal gases. We may 
summarize by saying that provided a suitable spectrograph is available the 
OH reversal method is likely to be simpler and more accurate than the Na reversal 
method for very hot flames. 
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ABSTRACT. A general process, by which light can be absorbed in insulators at frequencies 
less than the frequency corresponding to the energy gap, is considered. Using an Einstein 
model with a constant energy gap, calculation of such a process is made by means of the 
method of intermediate states. The results are tested for the case of cadmium sulphide, 
and agreement with experimental results is found to be good when allowance is made for 
mathematical approximations made in the calculations. 


Sil IN TROD WUCALLON 

WO previous attempts have been made (Moglich and Rompe 1940, 1942) 
| to explain the variation with temperature of the long-wave absorption edge 
of insulators. Neither method has given values, when tested, comparable 
with those obtained experimentally. In their first paper Moglich and Rompe 
proposed the method of multiple collisions, by means of which an electron 
interacts with several lattice waves and is lifted to the conduction band from the 
lowest band. Using the Debye model they made a calculation of the probability 
of such a transition. When this theory was applied to the shift of the long-wave 
absorption edge (Moglich, Riehl and Rompe 1940) a negligible shift, compared with 
experimental values, was obtained. In their second paper Moglich and Rompe 
suggested that the absorption edge would shift toward the red with increasing 
temperature, if allowance was made for the contraction of the energy gap between 
the first and second zones due to the expansion of the lattice. Calculation of this 
effect was made using the one-dimensional model of Kronig and Penney (1931), 
but the coefficient of linear thermal dilatation which was needed to produce the 
observed shift of 1A per degree K for the case of cadmium sulphide was 10-4, 
which compares unfavourably with the experimental value of 4 x 10~® obtained by 
Seiwert (1949). This theory is attractive, however, in that it may also be used to 
account similarly for the shift of the absorption edge to the blue with increasing 
pressure (Héhler 1949). In the particular case of cadmium sulphide, chosen 
because most experimental work has been done on it—by Seiwert (1949), Hohler 
(1949) and Kroger (1940) in particular—a simple consideration shows that, if the 
shift of 0-024 per atm. found by Hoéhler is completely accounted for by this 
process, it is impossible to obtain a thermal shift comparable with those obtained 
experimentally, namely 1A per degree K by Kroger and 1-44 per degree kK by 
Seiwert. For, consider a unit cube of cadmium sulphide, then a volume change 
of 1-5 x 10-* cm? due to the change in pressure of 1 atmosphere produces a shift of 
0-024. Thesame volume change, and hence the same change of energy gap under 
Méglich and Rompe’s theory, is produced by varying the temperature by 0:88°K 
when the absorption edge shifts by only 0-124. A further process must, 

therefore, be operating to produce the larger thermal shift. 
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In this paper we shall discuss possible transitions across the energy gap which 
are forbidden in the simple optical theory of solids. These transitions are per- 
mitted if the vibration of the lattice ions is taken into account. Considering the 
combined perturbation, we shall show that electrons from near the corners of the 
lowest zone can be raised to near the centre of the second zone, giving rise to an 
energy equation which permits light, of a lower frequency than allowed by the 
simple optical theory, to be absorbed. A linear shift with temperature of the 
absorption edge is also found. 


§2. PERTURBATION CALCULATION OF THE INTERACTION OF 
LATTICE VIBRATIONS AND ELECTROMAGNETIC RADIATION 
ON ELECTRONS 


In order to perform this perturbation calculation certain assumptions have 
to be made at the outset. It must be pointed out, however, that these approxi- 
mations are only introduced to make the calculation manageable and do not alter 
the basic physical ideas. Firstly, we shall assume that the lattice is simple cubic, 
with two electrons per atom, and that the zones do not overlap. ‘The first zone 
is then full and no conduction is possible. We assume also that the energy gap 
between the first two zones is independent of the temperature, thus ruling out in 
our model the process considered by Méglich and Rompe. For lattice vibrations 
we take an Einstein model. The wave function for the whole system is then the 
product of Bloch wave functions for the electrons and vibrational functions for the 
atoms. ‘Ihe quantum numbers are given by the wave vector k for the electrons, 
and n for the atoms, n having components in the x, y, z, directions. Both wave 
functions are normalized over an atomic cell. 

Expanding the perturbed wave function in terms of the initial wave functions, 
we solve for any one coefficient a(k,n,t) by the Dirac method of variation of 
constants. Fixing the initial state as kp, no, the following equation is obtained 
for any perturbation H,: 


., 0a 
— th = (Kemet) = (ke, me] A, | ko, ng) exp (El kong] — Elkeny])it/h, ee... (1) 


where (k;, n¢| H,| ko, no) is the matrix element between states ky, np and k,, n,;, and 
E[k,n] is the energy of the system in state k,n. In this problem the perturbation 
H, will contain two terms, caused by the light wave and the lattice vibrations 
respectively. It is given by 


H, =(eh/mci)A grad —R. grad V, 


where the vector potential of the light wave A =c€ [exp 27ivt + exp (—2zivt)]/47, € 
being the electric vector, v the frequency of the light, R the displacement of the atom 
from its equilibrium position and V, which is a function of the electronic coordinate 
r only, the potential of the unperturbed lattice. Inbothcases grad implies differenti- 
ation with respect to the electronic coordinates only. The displacement R must be 
taken sufficiently small for this formula to be true, so that the change in V is 
negligible outside the atomic cell concerned (Mott and Jones 1936). Without 
loss of generality the electromagnetic vector can be taken in the z direction, when 
A has the components (4,,0,0). Substituting from eqn. (2) into eqn. (1), it is 
easily seen that the matrix element for the transition as a whole does not exist, 
since the selection rules make one or both terms in eqn. (2) negligible or identically 
zero. ‘I'o solve the problem we must take a higher approximation in the solution 
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by Dirac’s method of variation of constants. This involves the use of intermediate 
states (Heitler 1944). 

Consider any intermediate state k,,n,, chosen so that the matrix elements for 
the transitions between it and the initial, and between it and the final states are 
non-zero, ‘Then, fixing kynp as the initial state, an equation similar to (1) is 
obtained with k,n,, replacingk,n, throughout. This equation will only involve one 
term of H,, the other giving a matrix element which is negligible or zero. Solving 
this differential equation for a[k,,n,,¢], the Dirac variation method is then re- 
peated with these values as initial conditions for a transition from k,n, to k,ng. 
This last calculation will involve the term in H, which was zero in the transition 
kon, to k,n,. The following equation for a[k,,n,,¢] is obtained: 


alk, ne, t= = (Ke Mal Ha] bo, mo)(Ker, ne Hi] kx, ns) 


k,n, (El kono] = E{k,n,]}) 
: {exe {(Elkono]—Elkene]+hy)it/K}—1 exp {(Elk.n,]— Elken¢] + hv)it/h} — i 
(El kono] — E[ kenge | + hv) (E(k,n,|-- E[ken,¢]+ Av) 


in which the time-dependent part of H, has been integrated between the limits 
0 and t, and we have neglected terms which lead to energy equations corresponding 
to the emission of light. Using the selection rules for the matrix elements only 
two distinct values of k, exist, and combining with each of these there are two 
values of n, for each of the three directions of vibration of the atom. The two 
values of k, are shown graphically in fig. 1, where E[k] is the energy of the state k 
which lies in the plane of energy discontinuity between the first and second zones 
and has its origin at the centre of the plane A where the energy is Fy. 


et 


A k 


Fig. 1. Optical transitions across the energy gap from an initial state k, to a final state kp through 
intermediate states k,’ and k,’. 


From the figure it is obvious that an effective energy gap hv is obtained which 
is considerably less than the true energy gap Avy, with corresponding shift of the 
absorption edge to lower frequency. 

The matrix element (k,”,n,|R grad V|ky, ny) divides into three similar terms 
on expanding the scalar product. Since V is symmetrical in the electronic 
coordinates, this may be written separating the lattice and electronic terms as 
3X, ng+i(k,,”|OV/dx| ko), where the x component of n, must change from ny to 
Ny +1 to satisfy the selection rules of the harmonic oscillator. ‘This reduces the 
number of terms in the summation of eqn. (3) from 12to4. Forming |a[k,, ng, ¢] |? 
we find that of the sixteen possible terms only four remain, the cross products of 
the terms in the sum disappearing because for large t they have no resonance values 
in the range of v, v<vp. We shall only consider one of these remaining terms in 
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detail, obtaining the results for the others by symmetry. Taking the term 
corresponding to the transition through the intermediate state k,”,% +1, and 
rearranging, we obtain the expression 

ks) 


(K sa") « (Ha ale 


Heine Ox 
¥ pee 1—cos(¢t/h) eee 
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6|X,,,. weal 


+O ° (+e — dere? 
in which ¢=£E[k,]—E£[k,”]+ W[mj)+1]—Wln], ¢=Elk,"]—Elke]+hv and 
Ay =(eh/mciv)(cl/27)"”, I being the intensity of the light wave which is defined as 
the energy crossing unit area per unit time. Wn] is the energy of a linear harmonic 
oscillator in the state n. 


§3. THE INTEGRATION OF EXPRESSION (4) 

The expression numbered (4) is integrated over all occupied initial and all 
possible final states. For insulators the lowest zone is completely full and all 
other zones completely empty initially. The lowest zone for a simple cubic 
lattice is a cube, while the second zone is the volume between it and the dodeca- 
hedron surrounding it. Using reduced wave vectors the second zone is reflected 
into the first zone (Mott and Jones 1936), and thus the region of integration for 
both initial and final states is over this cube. The integrand (4) is symmetrical 
in ky in any prism which is similar to one whose corners are given by (0,0, 0), 
(z/a,0, 0), (a/a, 0, z/a) and (z/a, 7/a, 7/a), where a is the lattice constant, and the 
origin of coordinates is at the centre of the cube with the axes passing through the 
centres of the faces. We shall consider only the integration of one term in (4) 
since the method of integration is the same for each. ‘Taking the third term we 

0 


have 
i= Deen ee wr od k,” 
= No, M+1 \, a | of J. ( ie 0x ag 
OV 2 1—cos[(e+ {)t/h] dk 5) 


(27)° 
zg (i. a ke) ellen te 


where Q is the volume of an atomic cell. The limits represent the prism whose 
coordinates are given above, the ky integration being a volume integral. To 
evaluate this integral the two matrix elements under the sign of the integration 
must be estimated. ‘Taking wave functions ,,=exp(k.r) sin {2mE,/fh?}?x, 


and },,.=exp(¢k.r) cos {2mE,/h?}"«, where k is measured from the point A,. 
we can evaluate the first matrix element. These wave functions are obtained from 
the nearly free electron wave functions (Wilson 1936) by neglecting a term which 
is small near the energy discontinuity. To this approximation the matrix 
element , is given by 

Ae k.") 


M?= (K: 9 Ox 


which is constant when k,;andk,” refer to the same k, otherwise it is identically zero. 

Using the Schrodinger equation, the matrix element (k,,”|@V//d«| ky)? reduces 
to |ky—k,”?| Vo—£?, where |V)—E| is the kinetic energy of the electron in its. 
lowest state at the surface of the polyhedron (Mott and Jones 1936). In order to 
simplify the integration of eqn. (5), we take an average value of |ky — k,” ?, and then 
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denote the matrix element by M,. Assuming that the average has been taken, 
but leaving the actual calculation of values until later, since the integration imposes 
conditions on ky and k,”, eqn. (5) becomes 


_ 28802 “A 2, [tel —cosl(e+ Oe/A] 
fo 1 ~cos[(e + ¢)i/] 
SOO ET, PMG! | * dk 2 { ake | dk, ig Toe es 


Taking the electronic energies of the intermediate and final states to be 
E{k,"]=(27+3?+ 2") and E[k,]=2,+«(£,¥2—x«)?+y?+2?, where Fp, is the 
lowest energy in the second zone, «=1+4E,/AE, and where x=k,h?/2m, 
y=k,h?/2m, and z=k,h?/2m. From this definition of E|k,”] and E[k,] we see 
that ¢ is a function of x only. Changing the variables from k,,, k., ko, k,, to (e+ 2), 
2, ky, x eqn. (6) becomes 


1440? 2m) 3/2 
mea (27/6 [Xn metal My ome {orl 
ee Sf > 1—cos[(e+ 2)t/h] d(e+ Z) 
dx | dky | dz | — eee 
ae of, a ler D{ater Qe” 2 
where a=Elk)] —«(£,1?—x)?— 22 +hv+hw— Ep, 
b= Elko] —«(£,1? — x)? —22?+hv+hw— Ep, 

and ha = W[n, +1] — Wino]. 


For sufficiently large ¢ the integrand can be regarded as a 6-function about 
(«+¢)=0, provided £=0 does not come in the range of integration. ¢=0 means 


0 LM N 3B 
FE {ko} 
Fig. 2. Region of integration of z and ky for fixed x (shown shaded). 


that v = v, at least, and since we are only interested in values of y <p this singularity 


does not occur. Similarly « cannot be zero in this integration if ¢ is not zero, 


since «+¢=0. Thus the limits may be taken as + o and — oo respectively 
provided a>0 and 6<0. Integrating eqn. (7) we have 


4 2m) 32 at 
I= ap Maem MEME EF 
a’ dx dz 
Pola’ k Sarr. “oe eT. NG Gee Poo Ee. ONO eoeevevee 
in Fo) ae | (Eik,] —a(E,i@— x)? Ey thy tho — 2)? °) 


where the limits of k, and z are determined by means of the original limits and the 
above inequalities. Putting a=0 and b =0, the inequalities define the region in a 
plane, x=constant, between the two full curves in fig. 2. The dotted lines are 
the conditions imposed by the limits of integration in eqn. (7). ‘The area over 


which the integration has to be taken is shown shaded. 
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Performing the integration over z in two parts corresponding to the ranges 
LM and MN of E/k,] respectively we find that eqn. (8) becomes 


_ 1440? (2m) 3? Pry eat) EA” dx 
= ay ea Xn, mors Ma Me's J C(x) 


7 yo a 7 
* | Rot | {sin ‘ DRE ) 5 i} ak | j 


The volume integration over k, is, for convenience, taken between surfaces of 
constant energy. If dS is an element of such a surface, the volume dk, becomes 
dk, =dEdS/|gradE|. Since the integrands are functions of E[k»] only, the 
integration over dS can be expressed in terms of the density of states N(E). 
Integrating over dS eqn. (9) becomes 


> 144Q (2m) 32 Nato ae ate Ey de B+ a 7 
= Gy ey Menkes | gL], NEM 


+ ae M(E) {sin a ae 
Be [Eko] —o( £7? — x)? — Ey thy + hw}? 4 J 


where B=a(£,1?—«x)?+E,—hv—hw. Since the range of integration of F is 
small and N(£) is a slowly varying function with no singularities, we may take N(F) 
constant over the integration. ‘The value of N(£) is adjusted to the average value 
of N(E[k,]) for the whole range of E[k,)]._ Performing the integration over E, 
eqn. (10) simplifies to 

1440 2m) 3/2 9S oe ear eee 

=o TEP Kanuri Pane me | eye (11) 
where C(x) =x —a(£,1? —x)?+hv— Ep. 

The other two terms in the expression (4), if treated similarly, give contri- 
butions for ¢=0 and «=0 respectively. Since we are only interested in hy less 
than AZ, the term for which £ =0 is omitted, and since the process we are consider- 
ing requires «+ ¢=0, we cannot have «=0 unless £=0, which implies v= Vo- 
Hence we can exclude the other term in (4) giving rise to the condition « =0. 

In evaluating integral (5) it will be noticed that it is immaterial whether we take 
«= E{k)]—E[k,”"]+hw or «=E[k)]—E[k,”|—hw. We now assume that the 
matrix elements satisfy the following conditions, (k,| 4,0/d«|k,,”) =(k,’ | A0/Ax| ko), 
(k.,” |AV//02| ke) = (k;V/Ax]k,’), and X,, =X, 4,1 approximately. Defining . 
the integral of a(k;, mn) overall initial and final states to be P(ko, k,)t, where Pis the 
probability of a transition from state k, to state k; in unit time, we find that 


576Q 2m) 312 or (Ha 2 
k,, kp) = =~ NM 222” 
P(Ko Ke) (27)3 (®){ zt [Xing m+1P My? Me j i Chee ee (12) 


The last integral is readily evaluated. 


§4. EVALUATION AND COMPARISON WITH EXPERIMENTAL VALUES 
OF THE ABSORPTION COEFFICIENT K 

Seiwert has determined the true absorption coefficient K for cadmium sulphide 
from measurements of the intensities of the incident and transmitted beams through 
a crystal of known thickness, reflection at the surface being taken into account. 

The relation between K in cm™ and the probability P of the absorption of a light 
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quantum of frequency v by an atom in time ¢ is given by K = NPhv/I, when N is 
the number of atoms per cm' in a crystal, and J is the intensity of the light at the 
point in the crystal to which P refers. 

In order to calculate P the matrix elements |X, ,,,,| and M, have to be esti- 
mated. |X,,,,,4,[° is evaluated in the usual manner by averaging the squares of 
the amplitudes of the oscillations of the atoms over all initial states. [X,,, ,,,,1[? is 
then given by é 2 THe ty 

mietil — qegip.....-. . ona 


where © is the mass of the atom, © the Einstein temperature, T the absolute 
temperature, and R Boltzmann’s constant. 

M,° was defined in §3 as the average value of |k)—k,”|? |V,— |? taken over 
all permissible values of ky and k,”, which we now know are those satisfying 
e+€=0. Since |ky—k,” |? =hy? +k," — 2hok,” cosa, where « is the angle between 
the vectors k, and k,”, the average over all possible values of « makes the third term 
vanish. ‘The average of k,2+k,” over all k, after substituting from «+ €=0 for 
Ry is (m/h?)[(14+«)E4/5+2(E_—hv—hw)]. Thus assuming a value of |V)—E| 
we determine ™,”. 

So far this calculation is quite general. In order to test its value, we shall 
compare the results obtained from it with the experimental values of Seiwert; 
a copy of his values are given in fig. 3(6)._ From his graphs we see that the total 
absorption at absolute zero would take place approximately at vy) = 5-906 x 101*¢/s, 
that is an energy gap of 2:43 ev. This limits the value of E, since we have taken 
the lattice to be simple cubic, because the highest energy in the lowest zone (3F',) 
must be less than EF, +/Avp) in order to make the substance insulating. We have 
taken FE, =1-12ev and FE, as 3-55ev. In addition, in an attempt to see if a Debye 
model would make any significant difference to the result, we allowed hw to take 
the value of the maximum of the Debye spectrum at the temperature considered. 
Taking |V,—£| to be 5ev we yive in the table the values of M,* at various 
frequencies and temperatures. 


Values of M,” x 10° at Various ‘Temperatures and Frequencies 


hv (ev) 
Di Dep 8 yy 
ee et 38°85 38-21 37°57 36-96 
T (° Kk) =654 : 
hw (ev) =0-103 38°43 37°83 37°21 36°61 
7) pede 38°23 37:63 37:03 36:41 


hw (ev) =0°135 


The error in taking the values of M,” to be the same as those obtained by 
replacing hw by —hw ine+¢=0 is negligible. ‘The values of K calculated from 
the values of the parameters given above are shown in fig. 3 (a). 

From these graphs the long wave absorption limit is seen to shift to the red by 
0-63 a per degree K. This compares favourably with the experimental value of 
1-44 per degree k when allowance is made for all the approximations included in 
the calculation. Comparison of the two graphs, however, shows a discrepancy in 
the rate with which the curves become asymptotic to the line v=vy.__ If this line 
were to shift to lower values of v with increasing temperature, agreement would be 
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considerably improved. This could be effected by a process such as that suggested 
by Méglich and Rompe in their second paper. We must conclude, however, 
that the process considered here is the most important since it gives good approxi- 
mations to both the shape and shift of the long-wave absorption edge. 


.; ray? ee 
x10" x1o'4 


Fig. 3 (a). Theoretical values of K. Fig. 3 (b). Experimental values of K. 
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Oriented Deposits on Crystalline Substrates 


By D. W. PASHLEY 
Physics Department, Imperial College, London 


Communicated by G. P. Thomson; MS. received 27th July 1951 


ABSTRACT. Layers of silver and thallium halides have been evaporated, in vacuo, on to 
cleavage surfaces of rocksalt, potassium bromide, magnesium oxide and mica, at various 
substrate temperatures. ‘These layers have been studied by means of electron diffraction. 
Examples of both high and low inter-atomic misfits between the substrate and the oriented 
overgrowth are found. In the case of silver chloride on potassium bromide, a chemical 
reaction occurred between the substrate and the deposit. Thallium chloride crystallized 
in an abnormal rocksalt type structure on potassium bromide, but had its normal caesium 
chloride type structure on the other substrates used. 


§1. INTRODUCTION 

URING the course of preparing specimens of various substances for 
certain electron diffraction studies, some results concerning crystal 
orientation on various crystalline substrates have been obtained. Since 
such results are of current theoretical interest (see van der Merwe 1949, Smollett 
and Blackman 1951), it is considered desirable that these data should be added 
to the existing information. ‘This paper therefore presents some results obtained 
from an electron diffraction study (reflection technique) of evaporated deposits 
of silver and thallium halides on freshly cleaved surfaces of single crystals of 

rocksalt, potassium bromide, magnesium oxide and mica. 


§2. SPECIMEN PREPARATION 

The specimens were prepared by evaporating the halide (zn vacuo) on to the 
various crystalline substrates, the halide being heated in a small tungsten spiral. 
The distance between the tungsten spiral and the substrate was about 12 cm, and 
the residual gas pressure lower than 10-*mmHg. ‘The substrate crystals were 
contained in holes in a brass block, during the deposition, and this brass block 
was heated by means of a coil placed below it. In this way the substrate 
temperature could be varied from 15°c to 250°c. ‘The thickness of the halide 
layers was usually of the order of 100 A, as determined from the amount of halide 
evaporated from the tungsten spiral. 
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In most cases the evaporated layers gave rise to diffraction patterns which 
consisted of one or more simple cross-grating patterns. These allowed the 
orientation of the halide to be determined. ‘The orientation of the substrate 
relative to the electron beam direction was given by its Kikuchi line patterns, 
which were still present, although rather diffuse. No Laue spots were obtained 
from the substrate, except when part of the surface was free from the deposit. 
The results of the interpretations are given in the table, and two of the more 
interesting of these interpretations are given below. 
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(1) Szlver Chloride on Potassium Bromide 


It was found, when silver chloride was deposited on to potassium bromide, 
that a chemical reaction between the deposit and the substrate occurred, so that 
silver bromide and potassium chloride were both formed. Fig. 1 (Plate) shows 
the pattern obtained from such a layer, which was prepared with the substrate 
at room temperature. ‘The Kikuchi lines arise from the substrate, and the 
analysis of the spot patterns is indicated in fig. 2. Parallel orientations of 
potassium bromide, potassium chloride, silver bromide, and silver chloride 
occur simultaneously. ‘The Laue spots from the substrate are observed in this 
case, because part of the surface was shielded from the silver chloride during the 
deposition. There is also parallel oriented metallic silver present; this formed 
during the electron bombardment of the silver halide (Pashley 1950). 

When the substrate was maintained at a higher temperature (150°c) during 
the deposition of the silver chloride, the diffraction patterns indicated only silver 
bromide, and no silver chloride or potassium chloride. The silver bromide 
pattern could then be identified both from its spacings and from the alternations 
in spot intensity which occurred due to structure factor. Such alternations are 


® Spots from KBr Substrate 


= “ 2 + + Spots from a parallel 
; drientation of AgCl 
° ° ° 

D r) © Spots from a parallel 

> a + jrientation & KCI 
e 2 Spots from a parallel 
bs Oar 8 8 om 2S paced real  AgBr 
e Spots from a parallel 

C ba yi 2 orientation of Ag 

> 


Fig. 2. The analysis of fig. 1. 


less marked with silver chloride. Under these conditions the silver bromide 
was either in the parallel orientation, or in a (111) orientation (see table). The 
two orientations occurred in very varying proportions from specimen to specimen, 
even when the two specimens were prepared simultaneously. 


(ii) Thallium Chloride Deposits 


When thallium chloride was deposited on rocksalt or mica, the normal bulk 
structure of thallium chloride (cubic, CsCl type) was observed. On potassium 
bromide, however, the thallium chloride was found to crystallize with a rocksalt 
type structure having a cube edge of 6:17+ 0-06 a. Fig. 3 shows the pattern from 
such a layer, with the electron beam set along the potassium bromide [001] 
azimuth, the strong cubic cross-grating pattern indicating that the new structure 
is in parallel orientation on the substrate. ‘This abnormal structure was usually 
accompanied by a small amount of thallium chloride with the normal structure, 
some of which was in the same orientation as observed on rocksalt substrates. 

Two further features are of interest. Firstly, the diffraction spots from the 
thallium chloride layers on rocksalt were often split up into pairs of subsidiary 
spots. This splitting was consistent with the presence of (111) type facets 

3-2 
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perpendicular to the surface, causing extensions of the reciprocal lattice points 
perpendicular to these facets. The identification of the facets cannot be ~ | 
considered conclusive, however, because the magnitude of the splitting did not 
allow sufficiently accurate measurements to be made. 

Secondly, the diffraction spots from the thallium chloride layers on mica 
were also split up into pairs of subsidiary spots, on two of the specimens examined. 
This splitting is illustrated in fig. 4, which shows the pattern obtained at the 
thallium chloride [211] azimuth. In this instance, the paired spots are due 
to a disorientation effect, and indicate that the thallium chloride (111) planes 
are rotated by approximately two degrees from being parallel to the mica (001) 
planes. Since the splitting occurred at all azimuthal settings of the specimens, it 
follows that the rotations occurred about all axes lying in the halide (111) plane. 
The [111] axis of any given thallium chloride crystallite thus lies on a cone of 
semi-apical angle 2° (see fig. 5), the distribution of [111] axes around this conical 


4 Perpendicular to Mica (001) plane 
| Alli] direction of a given TIC! crystallite 


Approx 2° 


Mica (001) Surface 


Fig. 5. The disorientation of thallium chloride on mica. 


surface being uniform. ‘The diffraction patterns from these specimens are 
obtained by considering each reciprocal lattice point to be replaced by a circular 
loop. No explanation of this effect has so far been found, but the effect might 
have some relation to the ‘spiral growth’ described by Frank (1949). The 
diffraction effects to be expected from crystals containing spiral dislocations have 
been discussed by Wilson (1948) and Frank (1948), and these effects bear some 
resemblance to the observed circular loops. This point will be discussed in 
more detail in a future publication. 

Fig. 4 also shows a diffuse streak pattern in the background. These streaks 
arise from continuous rods of intensity through all of the reciprocal lattice points, 
and along all of the [100] type directions. Similar diffuse patterns also occurred 
with the silver halides. ‘These effects will be considered in detail in a later report. 


§4. DISCUSSION 

‘The inter-atomic misfits between the substrate and overgrowth lattices is 
indicated in the table for all of the observed orientations. Although there are 
several examples of small misfits, there are also cases of fairly high misfits. In 
particular it is noted that silver chloride orients on mica, with a misfit of — pia sai 
Royer (1928) has grown many substances with a rocksalt type structure on ae 
by crystallization from solution. He found that substances for which the mish 
was less than about 15°, were oriented (with the same orientations as now 
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observed for silver chloride), but that for larger misfits no orientation occurred 
The fact that the silver chloride orients, despite the —24°/, misfit, therefore 
suggests that the method of growth of oriented deposits is an important factor 
to be taken into account by any theory of the process; an upper limit of percentage 
misfit observed for one method of growth does not necessarily apply to other 
methods of growing the same series of substances. Thus although the percentage 
misfit does appear to have some significance in connection with the appearance 
of oriented overgrowths, a small percentage misfit is certainly not an essential 
condition. The orientation of the silver halides on magnesium oxide presents an 
interesting case. A good numerical fit (misfit of approximately 6°) could be 
obtained if the orientation indicated in fig. 6 were to occur. This orientation would 
require ions of similar sign to be adjacent (e.g. at P, fig. 6), however, and this 
probably explains why it does not occur. 

When silver chloride was deposited on potassium bromide at 150° c, diffraction 
patterns due to silver bromide only were observed, so that all of the silver chloride 
had presumably reacted with the potassium bromide. In such cases no patterns 
due to potassium chloride were obtained (potassium chloride was only observed 
with depositions at the lower substrate temperatures). This might mean that the 
reaction occurred by movements of halide ions only, so that the potassium 
chloride was formed inside the potassium bromide substrate, possibly forming 
a mixed crystal with it. 


a ee AgCl [010] 
a ° 4 ° 2 < 
. fos 7 “Agel [001 
x 


@ Magnesium ions of magnesium 


‘N Z 
ee eee ey ee oxide (100) plahe 
we iMacs o Oxygen ions of magnesium 
aes zs ve RP, oxide (100) plane 
DR “™ > 
ee a - 7 3 > 
4 ee e Silver ions of silver 
) eo ° ‘eo ° chloride (100) plane 
vas (eS + Chlorine ions of silver 
4 ~*~ a \ chloride (100) plane 
4 
. ‘N 4 
¢ aN 


\ 
° » ° nis Mg0[010) 
y y « 
Lattices fitted 
here Mgo [001] 


Fig. 6. A non-observed orientation of AgCl on MgO. 


The final point worthy of note is the change of structure of thallium chloride 
from its normal caesium chloride type structure to a rocksalt type structure. 
Caesium chloride itself is known to switch to a rocksalt type structure above 
445°c, and the ammonium halides also undergo similar transitions. Recently 
Schultz (1950) has reported a case of caesium bromide (normally caesium 
chloride type structure) crystallizing with a rocksalt type structure when evaporated 
on to a lithium fluoride substrate. Such transitions must be influenced by the 
substrate, since the thallium chloride had its normal structure when grown 
on mica and rocksalt, but had the rocksalt type structure when deposited on 
potassium bromide under identical experimental conditions. ‘The fit of the 
abnormal thallium chloride structure on potassium bromide is only 6°, compared 
with 17°, for the normal structure. On mica and rocksalt substrates, however, the 
best possible fits which could occur for the abnormal structure are larger than 
those which actually occur for the normal structure. 
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On the Relations between the Photo-Elastic Properties and 
the Raman Effect in Crystals } 


By 0. THEIMER 
Department of Mathematical Physics, University of Edinburgh 


Communicated by M. Born; MS. received 24th May 1951 


ABSTRACT. General relations are derived between Pockel’s elasto-optical constants 
l ; 
and the polarizability derivatives ao, ,, ( a used in Born’s theory of the Raman effect 


in crystals. Numerical data are given for the case of diamond. Using these relations 
a formula is derived for the intensity of the Brillouin components of crystals in terms of 
the elasto-optical constants. This formula avoids the simplifications implied in an earlier 
expression due to Leontowitsch and Mandelstamm. 


§1. INTRODUCTION 
LL properties of the first-order Raman scattering are essentially determined 
ne l—I’ Meares 
by the derivatives «,,, , i) (p, o, w=1, 2, 3) of the polarizability 
tensor, which describe the changes of the polarizability of the atom (.) 


(/ denotes cell index, k denotes the different atoms in the elementary cell of a 


compound lattice), when any other atom i 5) is slightly displaced in the direction 
ve from the equilibrium position. c 

The numerical values of these characteristic constants are so far unknown; 
their calculation presents in most cases considerable difficulties.. There are 
three possibilities of determining these constants: theoretical calculation with 
the help of quantum mechanics, determination of the intensity ratio between 
Rayleigh and Raman scattering, and calculations based on photo-elastic 
measurements. 

The theoretical calculation is confined to very simple systems, because of 
mathematical difficulties. It appears that in no case has it so far been carried 


* Now at Physics Department, Technische Hochschule, Graz, Austria. 
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out; a possible exception is the hydrogen molecule, which, according to a 
remark of Burstein and Smith (1948) (without quotation), has been treated by 
Wang. 

The intensities of the Rayleigh and Raman scattering are simply related 
only in the case of gases. Here, however, the experimental difficulties of 
obtaining accurate values for the intensity ratio are considerable, because of 
the extremely weak intensity of the Raman scattering. The situation with regard 
to liquids, on the other hand, is confused from the theoretical point of view. In 
liquids, the scattering with unchanged or nearly unchanged frequency is 
a mixture of Rayleigh and Raman or Brillouin scattering, to which the polariz- 
ability «,, and the polarizability derivatives «,,,, contribute unknown 
shares. But even when this fact is not taken into account there remain other 
difficulties, such as the size of the region within which the scattered light is 
coherent, which are not known. 

- Similar difficulties arise in crystals. In an ideal crystal there is no Rayleigh 
scattering if the wavelength of the incident light is longer than twice the lattice 
constant.* All real Rayleigh scattering, which again in most cases can hardly 
be separated from the Brillouin components, results from lattice imperfections. 
For this kind of scattering no reliable theory exists. 

Most of these difficulties disappear if one tries to derive the polarizability 
derivatives from the elasto-optical constants p;,, which are closely connected 
with them. The only fundamental limitation of this method arises from the 
fact that the elasto-optical constants are essentially macroscopic constants and 
that their number is usually smaller than the number of independent 


RR’ 


Bes taeay which are atomistic quantities. Only if the symmetry of a 
crystal is high enough can some of the «,,_, (Gg at least be determined 


separately. ‘This is, for instance, the case with diamond. 

In this paper general relations between the photo-elastic constants and the 
polarizability derivatives will be given. The resulting formulae will be used 
to calculate the polarizability and the polarizability derivatives of diamond 
from photo-elastic measurements of Ramachandran (1950). 

Further, a formula for the intensity of the Brillouin components of crystals 
in terms of the elasto-optical constants will be derived. This formula avoids 
the numerous simplifications made in the derivation of the formula of 
Leontowitsch and Mandelstamm (1931) which have been discussed in a previous 
paper (Theimer 1951). 


§2. GENERAL THEORY OF THE DIELECTRIC PROPERTIES 
OF-CRYSTALS 


A well-known and remarkable feature of the atomistic theory of the dielectric 
properties of crystals is the existence of the two different formulae of 
Lorentz—Lorenz and Drude, which describe the relation between the 
polarizabilities of the single atoms and the macroscopic index of refraction. 
In the present paper a similar relation will be investigated, namely the relation 

ome ad oe if ; 
between the polarizability derivatives «,, ,, ( bh! of the single atoms and the 
* This statement is only approximately true, since there is always a Rayleigh scattering 


resulting from the mechanical and electrical inharmonicity of the crystals. ‘This second-order 
Rayleigh scattering can, however, be neglected in most crystals, compared with the first-order effects. 
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macroscopic changes of the dielectric properties, which are described by the 
elasto-optical constants p,,. It is, therefore, necessary to discuss shortly the 
problems connected with the existence of two different formulae for the index 
of refraction, in order to avoid ambiguity. 

For the purpose of this paper the difference between the two formulae can be 
considered to be purely formal, i.e. depending on the definition. In both cases 
the total polarization P of the crystal is represented as the sum of the electric 


moments p @ induced in the single atoms 


P=Ep(,). ee |) 


lk 


The p(,). however, are expressed in a different way: 


. @) as, (,) En h@seand ae (.) = (.) F pbk ee (2.2) 


Here the x (.) and the « (.) are the atomic susceptibilities* and polariz- 


abilities of the single atoms, E the average macroscopic electric field in the 
crystal and F the effective microscopic field polarizing the single atoms. 
Definition (a) corresponds to Drude’s formula, (6) to the formula of 
Lorentz—Lorenz. 

This formal aspect of the problem is, of course, of little importance compared 
with the numerous questions of real physical interest which are connected 
with the determination of F. But it is just this formal aspect which is important 
for the present paper, as it allows the unambiguous definition of the physical 
entities used in the following sections. Since in all existing papers on the 


Raman effect in crystals the dipole moments p ‘ induced in the single atoms. 


are represented as functions of E and not of F, the definition (2.2 a) will be used 
which justifies the application of equation (3.7). It is, however, convenient to 


change the notation and to write x ( for the atomic susceptibility, since this 


k 
is the common usage in all papers on the Raman effect. 


§3. THE RELATIONS BETWEEN THE ELASTO-OPTICAL 
CONSTANTS pj; AND THE POLARIZABILITY 


DERIVATIVES «go, 1 (a 


Let « be the reciprocal tensor of the dielectric constant, which determines 
the index ellipsoid and S an homogeneous strain defined by + Soo =OU,/OP 9, 
Soo =0u,/O7, + Ou,/Or,. where r=position vector, u=displacement vector of a 
lattice point. p. o=1, 253. 

* The macroscopic susceptibility « is defined as the polarization per unit field and unit volume. 
The atomic susceptibility « h)? introduced here, is however not related to the unit volume. 


The notation has been chosen to indicate that the summation of z( ) over all atoms in 1 cm® gives 
immediately the susceptibility. k 

_ Usually Soo is defined as Soo=}(0ug/0rc+ Ouc/ Ore). Here, according to Pockel, the factor a 
is omitted; thus terms with mixed indices can then be treated in the same manner as terms with 
equal indices. 
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The homogeneous strain S produces a change in ¢!, which in first 
approximation may be considered as a linear function of the deformation: 


Aex, <= '=E Po urS fara weil ee (3.1) 


The Pyo,,» are Pockel’s elasto-optical constants. 
For such small changes in «+ and the dielectric tensor «, we obtain, by 
differentiating the identity 


% 9 iy ee Lae < Be cere (3:2) 
the first-order relation: 
Dike, ee Fen Nene See Oe qn Wenn, fee (3.3) 
1G 
For cubic crystals for which 
ew = bof au i =6,,€71 Seater (3.4) 
the above relation reduces to 
x NE 9 
ee erie ae i ee eee etek. (2.5) 


The corresponding changes of susceptibility Ax,, can be found in the 
following way: put (D=dielectric displacement vector) 


D2=<E=E+47P)and® P=cEw 7 st: (3.6) 

then Kee Lae ae ek ek Meera (37) 

Hence Bape re gAt = Lhe yO uss - 9,” eee eee (3.8) 
py 

where Ce MEAT) Pig ye 1 PR eee (3.9) 


: ; l : 
We now derive expressions for the «,, ,, 1n terms of the «,,_, ee which 
are used in the theory of the Raman effect. 


Consider first arbitrary displacements u(;) of the atoms al and denote 


with Aa, the changes of the susceptibility of the atom (.) , caused by 


l 
i ; 
the displacements u (.) of the other atoms. ‘Then 


l I-Il\ (l 
=m ae 
oe @ 7 foes i ( kk’ ) U,, i (3 10) 


The total change of the susceptibility of the whole crystal is 


l l-I’ I 
= =» gb oan d : 
jae edt (.) : 3 3 tea u( hk’ )u, ) (3511) 


In the case of a homogeneous strain S,,,=0u,/0r,+0u,/Or, the displacements 
are linear functions of S,,. ‘These functions, however, are different in the 
classical continua theory and in the correct atomistic theory of Born. In the 


classical continua theory they are 


i OU, Ny 
46(y)= Rae apt tak cae (3.42) 


In lattice dynamics (3.12) is a zero order approximation which holds 
only for the centres of gravity of the single elementary cells. ‘The actual 
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i @ of the different atoms in the elementary cell of a — 
compound lattice differ from (3.12) by a term independent of J’ of the form ~ 
le Ou 

(1) = ) coped ie “tea. Ameria 3.43 

10 (5) = BE al) Fe G.13) # 

In the general case I’,,,(k’) is a somewhat complicated function of the 
dynamical matrix of the lattice, which has been given by Begbie and Born (1947).* 


displacements u 


; L oe . 
The total displacement u( connected with a homogeneous strain is from ~ 


(3.12) and (3.13): : 
/ / rs) f iff f 0 . 
U, es =4 (a + u,0(k') == E r, (.) + LD A(R’) = |: re? (3.14) 
‘ y v q 


The susceptibility changes of the single atoms k depend, of course, mainly 


ay which are not too distant from 
them, i.e. the first, second, third... neighbours. To express this (3.11) applied 
to a homogeneous strain has to be written in a different form. Write 


(ji a5 
where u ee a) = u(/)—u(/')+u(k)—u(R’), i (3.16) 
Then (3.11) becomes 


IONE ee il! 
Mtge = EB aaa gy l= G)- ee mE Sn (3.17) 


Put /—J’=/ and use the relation X «,, , Gv) = which follows from the 
Ik’ 
invariance of «,, against rigid translations of the whole lattice: then 
(N=number of elementary cells per cm*) 


l l 
=-NiX 2X <a eer : 
Baye= NEE E tangy) (ay) ea 
Introduce (3.14) in (3.18), then 


a I Ou, l , du, 
IPE bg SUS (ww) Se aes +E Dan (RR) 
v 7] n 


: l 
on the displacements of the other atoms ( 


- ona (3.19) 


uv Kk Uk’ 
Rearrange with respect to the du,/dr,, then 
Ou 
Ad oo ae Xoo be a ee ee ee So. A (3.20) 


Hh a l l x l ns 
wheres a7, = sabe 2 eae G3) ee (5) a thea, 9 6 DT w(RR )| ‘ 


(0u,,/O7,)%oo, ny Tepresents the changes of the susceptibility %), Of the crystal 
connected with a shear du,/dr, of the angle between the p- and p-axes of the 
undeformed crystal. The same shear and the same susceptibility change 
results from a homogeneous deformation du,/Or,, if Ou,/Or,=Ou,/Or,. Hence 


(G4,,/O9, Joc ge, jig = (OUsIOT,\ om, «nih en ee (3.21) 


or Oh og. ivi Sgas use are ule ee (3.22) 

* (3.13) can be derived from eqns. (2.5) and (2.22) of Begbie and Born if one replaces gy b 
(1/ug) Cug/er,. This substitution is allowed for long elastic waves as can be seen from differentiate 
Ug= Ug exp (7q.. r). 
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With the help of the invariance condition (3.22) and taking into account 
that the factor $ has been omitted in the symmetrized strain components 


Si, = Ou,/Or, + 0u,/Or,,, (3.20) reads 


GN Sa ae eet er (3.23) 
ao 
A comparison between (3.23) and (3.8), (3.9) yields 
l l l : é? 
> = Eas Go) re (es) + 7 Cpe a (2 DP RR )| = Gr Peo ree (3.24) 
This is the general relation between the susceptibility derivatives «,,. , A 


and the elasto-optical constants p,, ,, in cubic crystals. 


§4. NUMERICAL VALUES OF THE (go, x (ae IN DIAMOND. 


The photo-elastic constants of diamond have been determined by 
Ramachandran (1950). There are only three independent p,,,,, in cubic 
crystals. ‘They are given (in the notation of Pockel) in table 1. 


Table 1 
Pu=Po0, 00=DPoo, co= Pup, yp= — 9°31 
Pi2=Po0, co=Po0, yp=Poo, eg=DPoo, yu =Pup, e¢=Pyp, oo = 0°09 
Pas=P oo, 00=Pop, ou Pug, wo = Poe, 00 -=Pyo, po =P oy, ox = — 9°12 
All the other foo, ny are zero 


The great number of «,, , fe) is reduced by the particularly high 


symmetry of diamond to four independent constants, if the simplifying 
assumption is made that the susceptibility of the atom k depends only on the 
displacements of the first neighbours and of the atom 2 itself (J=6, 1, 2, 3, 4). 
‘This assumption is quite admissible in atomic lattices with homopolar binding. 


SN RR 
derived in a previous paper (Theimer 1951) and are presented in table 2. 


The symmetry relations between the « ( 5 ) of diamond have been 


Table 2 

‘oa, H kk’ 0 1 2 B ; 4 Xoo, 7 kk’ 0 1 2} 3 4 

11,1 0 A —-A -A- A = %9,1=O93,3=O%g1,1= Maa, 0 FE —-E E —-E 
O99,2 OF A A A —A O39 31,1 = %g2,2= 13,1 One Fee He Ee 
33,3 0A A —A —-A_  %1,3=042,2=013,3=Oa1,2 Oe — fn Ee E 
O41,2—= 33,0 0 B -B B —B O19. 3—=MXa1,3 = 7:99) JD) D D D 
22,3 M113 0B BB —B  o43,1=%,1 TAD 1D Die ard) 10, 
33,1 M91 0 B —B —B B X31,9=%13 2 —4D D D D D 


Reto ces ) (I denotes the neighbours of the atom A’) are found with 


the help of the relation «,, , GS) == Koon (ee) . 


The values of 7, ie for diamond have been given by Smith (1948); 


they are contained in table 3. 
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The values of [,,,, are (Theimer 1951) 
D-P_(k 0-66 a/2 for pAvAn, 
k’ ws ip <= “WY \ = 
LAE eh) sige 0 for all other cases. 


Using tables 1, 2, 3 and eqns. (3.24) and (4.1) the susceptibility 
derivatives A, B, D, E are (R=e?/l6mNa=5-88x10-%, ¢=n?=6-81, 
n=2-61 =index of refraction) 

A=p,,R= —1:82x10-% cm’, B=p,,.R=0:53 x 10-% cm?, 
(E+0:66D) =pyR = —0-71x10-%cm*. —........ (4.2) 


For the susceptibility «(k) of the carbon atom in diamond one gets from. 
the definition (3.7): 


1 e—1 
=> = 2Natk) = 14.01 9 i ee 4.3 
i © “(f) 4a Ce) 
ag ue on ay ee 
a(R) = Sa N = 2, 62 x 10 CIS. ~— sy Vet © So pees iecraete (4.4) 


These values for A, B (H+0-66D) and «(k) give calculations of the relative 
intensities of the Brillouin components of diamond in good agreement with 
experiment (Theimer 1951). The numerical values given in (4.2) and (4. ui are: 
valid for radiation with wavelength A=2536 a. 


Table 3 


DE 
aK 
i) 
WwW 
aN 
pe 
eI 
Nl 
wl 
+| 


xy La ka a $a x,  —4ta 4a ta  —ha 
Xe La ta La ha X,  —ta La da da 
x3 a ta —ta —ta x3 —ta —tha ta 4a 


a=lattice constant=1:78 x 10-8 cm. 


§5. THE BRILVOUIN BEFECT IN CRYSTALS 


The intensity of the Brillouin components of crystals can be calculated from 
the following formula aa LO5 Ds 


ee. 


Fo) ian oO se Oyen wo(® at re eee Shen (5.1) 


e 
where A, A* are the complex amplitudes of the incident light wave, s=unit 


vector normal to the direction of observation, v Q is the frequency of the 
e€ 
lattice normal vibration j, with wave vector @, which fulfils the condition 


Q=-Q’—Q’. (Q’, Q” are the wave vectors of the incident and the scattered. 
light.) 


2) or { ak l 
aL =I Sak 
ns (, eas Sas le (2 Sait (3 


+ Bes fy Pak, | AG ) RARER Ls. (5.2) 


i) 
is the displacement of the atom & if the lattice performs the 


Here’ 1, (i | rs) 


vibration 


The index 7,=1, 2, 3 denotes three elastic vibrations, which 
e 
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can roughly be interpreted as one longitudinal and two transverse waves. 
Eqn. (5.2) can be expressed in terms of the photo-elastic constants. 


Multiply (3.24) by 7Q,u, (°) and sum over all yp, v, then 


l l l O 
yyy / 
ad _ i « Koa, “ Gs ) Os. Gay, ni 2ot9o, a Pru lkR 10 Ju, (i | 4) 


le” 
= FE Peo Qa (: 5) is a er (5.3) 


(5.3) is identical with the second term of (5.2). The first term of (5.2) is, 
according to (4.3), 


IND Zttgo() Qnty (2) =i 9 tee (A): ey (5.4) 
(5.3) and (5.4) Satine give 


no (2) 2 Bflen-D0t PDPamOrheg (2). 2-6) 
Q 


There are still the eigenfrequencies v c 


Je 
which can be replaced by macroscopic constants. In the region alQ|<1 


) in the intensity formula (5.1), 


(a=lattice constant) there is a velocity of sound v ) for each elastic branch /,, 


which does not depend on |Q| but only on the direction Q@/|O|=s; hence the 


frequency is 27v ice = o(; ) [QO]. 


The velocity of sound v 4). is completely determined by the elastic 


e 
constants c;, and the density p of the crystal. Hence the intensities of the 
Brillouin components can be represented as functions of the macroscopic 
constants ¢;,, p,, and». The simple form of (5.5) only holds in the case of cubic 


crystals. For crystals with lower symmetry the relation between the «,, Gs ) 


é 
and the p,,. ,,, 1s, according to the general formula (3.4), much more complicated. 


Further, the macroscopic formula does not allow the calculation of the intensity 
ratio between the Brillouin components and the optical Raman lines. This 
ratio can be obtained only with the help of (5.1) and (5.2) and the corresponding 
formula for the Raman effect with optical vibrations (Theimer 1951). 
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ABSTRACT. A theoretical investigation is made of the amount by which a dissolved 
gas should reduce the tensile strength of a liquid. The tensile strength is assumed propor- 
tional to the surface energy of a vapour nucleus, which is further assumed to have the form 
of a product of a function of its radius and a function of the dissolved gas concentration. 
The surface tension of a one-component system is expressed in terms of the van der Waals 
constants (developing Fowler’s theory) and the treatment is extended to a two-component 
system on the basis of Wall and Stent’s theory of van der Waals binary mixtures. The 
conclusion is reached that the reduction in the tensile strength of water caused by saturating 
it with air at a pressure of one atmosphere is less than $%. 


§1. INTRODUCTION 

HE present theories of the tensile strength of liquids all appear to give 
| values very much higher than experiment. Most of the experiments, 
such as those of the Berthelot tube type, give values for the tensile strength 
of water less than 60 atmospheres (Temperley 1947). ‘That this discrepancy 
is, at any rate partially, due to the experimental conditions is suggested by the 
recent experimental work of 'Trevena (1950), which gives values for water of up 
to 300 atmospheres. Similar values were also obtained by Briggs (1950), 
Ursprung (1915) and Renner (1915). The lowest theoretical value, due to 

Temperley (1947), is about 500 atmospheres. 

Many authors attempt to explain this discrepancy between theory and 
experiment by arguing that the ability of a liquid to withstand tension is reduced 
by the presence in it of dissolved gases; others believe that small bubbles of gas, 
in the liquid or in Griffiths cracks in the containing walls, are mainly responsible 
for the lowering of the tensile strength. ‘The present theory attempts to decide 
between these possibilities by estimating the amount by which the tensile strength 
of a liquid is reduced by dissolved gas. 

Temperley’s theory is based on the use of a van der Waals equation of state, 

(p+al?)\V—b)SaRT —. (> Eee (1) 
For suitable temperatures (less than some value 7\,) it is possible to obtain 
metastable regions of negative pressure, i.e. tension. 

On the other hand, theories such as that of Déring (1939), Firth (1941), 
Frenkel (1946) and Fisher (1948), which attempt to calculate the rate of growth 
of a gas bubble, predict values for the tensile strength of water greater than 
1000 atmospheres. Consideration of the stability of a microscopic bubble 
clearly necessitates a knowledge of the surface tension of the liquid. It is usually 
assumed that the surface tension at the boundary of such a bubble is the same as 
that over a plane surface, but, as the bubble must be of the order of molecular 
size (if its probable time of formation is to be less than 1 sec, say), this assumption 
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appears unwarranted. Since the energy needed to create a small bubble and the 
van der Waals constant a are both measurements of the intermolecular attraction, 
we shall assume a value for the effective surface tension of a small bubble that 
will give a predicted tensile strength in agreement with that given by Temperley’s 
theory. 

We will write o = o(r), the surface tension of a bubble of radius 7 and o° = o(00), 
the surface tension of a plane surface. Then for the bubble to exist in (unstable) 


equilibrium, RIG ks ER RE Rah, (2) 
where p and 9’ are the pressures in the bubble and liquid respectively. If the 
tension t=—p’ is greater than the value from (2), the bubble will expand 
indefinitely, and the liquid will rupture, while if it is less than —p’, the bubble 
will collapse. Hence the tensile strength is t=20/r (neglecting p which is small 
compared with 2). 

If o,, o,° are the corresponding quantities for the pure liquid, we assume 


Te OT ee we a SE ive ree (3) 


We shall also denote the tensile strength of the pure liquid by ¢,. 7, in the ‘ hole’ 
theory, is the radius of the largest bubble which can be formed by fluctuations, 
in a reasonable time, and is of the order of d, the intermolecular spacing. We 
assume 7 unaltered by the dissolved gas. ‘Then t, =20,/7, and therefore 


A GET yee ea ras ae (4) 


Hence the problem of finding the effect of dissolved gas on the tensile strength 
reduces, on our assumptions, to that of finding the effect of the solute-on the 
surface tension. ‘The assumptions seem justified if the arrangement of molecules 
in the liquid is not greatly altered by the dissolved gas, which is probably correct 
until the concentration reaches several per cent. We now proceed to discuss 
the surface tension in some detail. 


§2. SURFACE TENSION AND THE PARACHOR 
For a liquid in equilibrium with its own vapour the surface tension o is given by 
isle pee ieee. pe Ve oe,’ (5) 
where p,, py are the densities of the liquid and vapour respectively, M* the 
(chemical) molecular weight and P the so-called parachor. ‘This equation is 


due to Sugden (1930) and is discussed further in the Appendix. Let us write 
ol4—7, (P/M*).p,.=7,, (P/M*),p,=7,. Then in virtue of eqn. (5) we have 


rod re Mel a Fee © Ano (6) 


Thus if 7, and 7, can be determined then 7, and hence o, is known. 

The determination of +, and 7, for the liquid and vapour phase respectively 
demands a knowledge of the density, molecular weight and the parachor for that 
particular phase. The first two quantities can be readily found but the deter- 
mination of the parachor is not so easy. Using the results of a paper by Fowler 
(1937) it has, however, been found possible to express the parachor of a substance 
in terms of known quantities if one assumes that the substance obeys a van der 
Waals equation of state; this extension of Fowler’s theory is given inthe Appendix, 
where it is shown that the parachor P is given by 


PES ol Ee ee eR (7) 
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where w =[3(3«’ + 1)/647.N49]14, «’ =1-91, N=Avogadro’s number and a,6 are 
the van der Waals constants for the substance. 

Thus, for a van der Waals liquid and vapour, it should be possible to determine 
7, and 7, respectively and hence o. We shall now proceed to do this. 


§3. THE SURFACE TENSION OF A VAN DER WAALS 
TWO-COMPONENT SYSTEM 


We have seen that if the equation of state obeyed by a liquid is that of van der 
Waals then we can determine the parachor and hence the surface tension o in 
terms of the van der Waals constants. 

Let us therefore consider a gas whose van der Waals constants are a, }, 
which is dissolved in a liquid whose van der Waals constants are a, b,; we shall 
treat this mixture as a van der Waals binary mixture. Such mixtures have been 
discussed by several workers. 

If x, and x, are the molecular concentrations of the two components, then 
van der Waals assumed that the constants a, 6 for the mixture are given by 

A= X17, + 2% X%2Ayo + X90, b= 4,654 9505, nome (8) 
where a. is a constant depending on the inter-species interaction. In our case 
component | is the gas and component 2 the liquid. Let x be the concentration 
of gas in the condensed phase. Then x,=x, x,=(1—.~), and the two constants 
for the mixture in the condensed phase are given by 

A, = Ay + 2x( Ay, — a2), b= by + x(6b35), 2) Gees (9) 
where we have neglected x?, which is justifiable for concentrations of gas up to 
10%, 

Let y be the concentration of gas in the vapour phase. ‘Then, in a similar 

way, we find that the constants for the mixture in the vapour phase are given by 

Ay = Ay + 2V(Ay2 — Ap); 6 = 6, (0 —0.). ees (10) 
where y’ has been neglected on the ground that y” is small when x? is small. We 
justify this below. 

We must now proceed to determine 7, and then 7,. 


(1) The Condensed Phase 

In this phase, since R7<a,b,, then for one mole the volume V =), 
approximately. (This amounts to taking the point C instead of A in fig. 1.) 

The units of V are cm? mole; hence if p, is expressed in g/cm? we have 
p.= M*/b... Now for the condensed phase, the parachor P, is wa,¥/4b,72. Thus 
7 = Pb, =og"*bh = oF 

T. = Wg 4b. F241 + ¥(Ay9/2a, —56,/12b,—1/12)} ...... (11) 

in virtue of eqns. (9). ‘Thus 7, has been expressed in terms of the van der Waals 
constants of the two components of our mixture. 


(11) The Vapour Phase 


Wall and Stent (1949) give the following equations for the partial pressures 
of a two-component van der Waals mixture 


} d Roane 1: (ab 2X 1A, + 2%X.a 

oper. | dp 2) 9 J eas 1 $40, XQ, 242 

SPi gx og (a,/b, yar { Ae b ye RT b2 b.. 5) 
ui tests (2) 


=log x, +loe(a,/b2) 4.2 2a tMdis 2) | 1 facb, — 2itga,+ 2x dys 
log po § gr g ( /O.") +2 + eee ee 
7 j il ae b. th Gee b. y 
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where x,, x) are the concentrations in the condensed phase. From eqn. (12) we 
have, after some reduction and substituting for x,, x, and a,, b., 


log (p,/pP2) =log x + {2 (Se — 3) + ERT ( 2a. + ay) | + O(x) 
=logx+ {a+ B/T + a,/b,RT} + O(x) (Say ieee eye co. (13) 
If we put C=«+8/T+a,/b,RT then we have 
Pies =—Serpuc Ole) a= | es (14) 


Since y is the molecular concentration in the vapour phase, then, provided the 
vapour pressure is small, we may apply Dalton’s law of partial pressures, i.e. 


Pi = PilP2 (15) 
Pr Pe LEE Di Payee SY aie tee te Sa 
Thus, in virtue of equation (14) we have y=xe* to the first order in x. Hence y 
is approximately proportional to x; in other words, if the concentration of the 
gas component in the condensed phase is small, so also is its concentration in the 
vapour phase. 


vies 


Fig. 1. Van der Waals isotherms. Provided Fig. 2. Molecular interaction potentials of 


T<Ty, the liquid exhibits a _ tensile the form E(r)=—p/r°+ A/r”. Dotted 
strength t. The broken line AB is the curve n=12 full curve n> «<. The 
equilibrium curve. latter form has been used in the text. 


From eqn. (9) we have (a,/b,”) =(dy/bs”)(1+ax). Also, from eqn. (12) we have 
a 


- a; 5 
pra 55 {exp («+ 7) fe + O(«)}, ?.=(1 —x)F3 {1 _ iRT —ax+ O(w*)| i 


the definitions of « and f being given in eqn. (13). ‘Thus the total pressure p in 
the vapour phase is given by 


P=P1 + Pa = (ay/b,?){x exp («+ B/T) +(1—x) exp(— a/b, RT)} 
neglecting O(x?). Since a,/b,.<RT, exp(—4@,/b,RT) may be neglected in 
comparison with exp(«+f/T). ‘Thus 


p= 5% {exp (—a,/b,RT)}{1 +.vexp(%+ B/T+ a/boRT)}.....-. (16) 
2 


As we are at the point B (fig. 1), we may apply the perfect gas equation with 
sufficient accuracy. We then have 1/V =p/RT or 
oe a a 2 TERT xP (—@/5RT)}{1 +wxexp («+ B/T +a/b,RT)}. 
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Then 7, is given by ered NAb aga n ene Sad RG eee (17) | 
Now | 
Ay2—4 7 b,—6 | 
4p 7/12 — 9 4p 7/12 aoe, LPO ; 
atthe ay"!4b, { 1+y ( oe +i e )} from eqns. (10) 


= ay 47211 + 0 ay9/2dy + 7b,/12B5 — 13/12) exp (a + B/T + ay/b,RT)}, 


since y=xe°. So, after some reduction, we have, on substituting in eqn. (17) 
arag th 2 —a? Gige 6 iby ol B Ay 
7 RTs PERT I. 1 ag 12byi 1) ee ee 


Thus 7, is determined in terms of the van der Waals constants. Hence, from 
eqns. (11) and (18), we have 


n=T.— a, 
a 56 1 a —a 
= 1/47, —5/12 ey ed et) Tees eat 
mee aUa. {14x (3H 12b, 3) a has 6 Re) 
dia os GAUe 1 
~ gape +h 7 (5 ot yp a) 


or, neglecting unity compared with exp(«+/T), 
Seep 1s pea 
Tt =wa,1!4, {1 RTO? (5 ee) 
a ae BN: pair eee 
5,RT exp («+ 7) | (x + (2b, <h1D Re ea (19) 


Thus the surface tension o° is given by 
4a, 
b,RT 


~ pRiplesp (@+8/Ty (8 +H -3)} +++4+(20) 


to about 5%, since we have neglected the square of exp (—4,/b,RT). 


ohm ash {1 ma exp (—4,/b,RT) 


§4. REDUCTION IN TENSILE STRENGTH OF A LIQUID OWING TO 
DISSOLVED GAS 


We can now find the effect of dissolved gas on the surface tension, and hence 
on the tensile strength, of aliquid. For a pure liquid x=0, and so from eqn. (20), 


we have o,,° =arta,b,53 {1 _ (4a,/b,RT) exp ( _ a,/b,RT)}. Sods (21) 


Remembering that t/t, =0°/o,° we have the following relation between the 
tensile strength of the liquid with dissolved gas and that of the pure liquid 


2a): 
t/ty =1—(wa,/b,.RT)[exp («+ B/T)] 3 | (22) 
a, 36, “3 
For small concentrations of gas we have already seen that 
Py = (%42/b2*) exp (% + B/T) 
At ibs (2a 70). 901 
and so we have fag ~ in (= cae mew ee (23) 


This equation gives the relative decrease At/t,, in the value of the tensile strength 
which is to be expected due to the presence of the dissolved gas. 
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§5. NUMERICAL VALUES 

Let us now take as an example the case of water containing dissolved air. If we 
assume that aj) = 1/(@,a,) we can compute the weakening caused by each of the 
component gases of air by means of eqn. (23). The data are listed in table 1. 
We then find, if the water is saturated with air at one atmosphere pressure, 
At/t,.=7 x 10~* for oxygen, Aé/t, =31 x 10-4 for nitrogen, Az/t, =0-3 x 10-4 for 
argon. (p,=0-8 atmosphere for nitrogen, 0-2 atmosphere for oxygen and 
0-01 atmosphere for argon.) 


Table" 
Substance H,O Os N. A 
Critical density (g/cm*) 0-329 0-430 0-311 0-531 
Critical pressure (atmospheres) 218 49-7 53-5 48-0 
a (atm. lit?) 0-490 0-0203 0:0202 0-0200 
6 (litres) 0:0274 0:0372 0:0451 0:0376 


Thus we see that the maximum weakening caused by air dissolved at a pressure 
of one atmosphere is about }°% and that the amount of weakening is independent 
of the solubility of the gas, provided only that this is small enough for the 
approximations to apply. The smallness of the decrease in tensile strength is 
rather unexpected. It should be emphasized that this theory gives the order 
of magnitude of the effect only, but it is unlikely to be changed significantly by 
the use of different van der Waals constants. 


§6. DISCUSSION 

The conclusion that the presence of dissolved gas only causes a negligible 
diminution in the intrinsic strength of liquids seems surprising, yet it does not 
seem to be in disagreement with the experimental evidence. For example, the 
high values obtained by Ursprung (1915) and Renner (1915) were for the tensile 
strength of cell sap, which cannot possibly have been air-free, and was, in all . 
probability, nearly saturated. Again, Scott et al. (1948) introduced a very 
elaborate technique for preparing Berthelot tubes containing water with a very 
small gas content, yet they did not find any significant increase in the observed 
tensile strength over that found, for example, by Meyer (1911), Worthington 
(1892) or Temperley (1947) in glass apparatus. ‘Thus, in accordance with the 
speculations of Vincent (1943), the discrepancy between theoretical and. 
experimental values of t appears not to be explicable in terms of dissolved gas, 
but must presumably be sought in the small bubbles of air trapped in Griffiths 
cracks in the wall of the container. 

On the other hand, there does seem to be definite evidence that the observed 
tensile strength, as measured by ultrasonic methods, does increase as the gas 
is removed (see, for example, Blake 1949). ‘This increase may be, at any rate 
partially, due to the removal of undissolved gas nuclei, and further the mechanism 
of rupture may involve the process termed by Blake ‘rectified diffusion’, that is 
a progressive diffusion of gas into microscopic bubbles in forced oscillation in 
the ultrasonic field. Such a mechanism is absent when the stress is purely 


hydrostatic. 
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APPENDIX 


THE SURFACE TENSION OF A DIETERICI LIQUID 


The assumptions made in deriving either van der Waals’ 


(ptal/V*)V—b)=RT -— sane (1) 
or Dieterici’s equation of state 
p(V—b)=RTexp(—a/RTV)  —_.... wae (A 1) 


are of rigid, perfectly elastic molecules, which are subject to a short-range 
attractive interaction. These are just the assumptions on which any simple 
qualitative explanation of surface tension is based, and it should therefore be 
possible to express the surface tension of a liquid which obeys such an equation 
of state in terms of the constants a and b of the above equations. The two 
equations are identical, provided we neglect O(1/V”). 

Macleod (1923) found empirically that o 0 (pjiq—pyap)* over a wide range of 
temperatures, where o is the surface tension, and the p’s are the densities in ~ 
g/cm*. Sugden (1930) defined the parachor P by the equation 


i, = M*a"4 (pig — Paap) 6 tee ens (A Z 


where M* is the (chemical) molecular weight of the substance. He suggested 
that the parachor might be a useful parameter for studying the structure of 
complex melecules. 

Fowler (1937), in a statistical theory of the surface tension of liquids, showed 
that Macleod’s formula should hold near the critical point. In this region, 
eqn. (A1) rather than (1) should hold, and we therefore follow Fowler in treating 
a Dieterici rather than a van der Waals liquid. But having found the surface 
tension in terms of the constants a and b, we revert to van der Waals’ equation 
in the text of the paper, where the conditions are far from critical. (A Dieterici 
liquid has no tensile strength, on Temperley’s theory.) Fowler derives the 
following expression for the parachor: 


aN? M* (° 
Pear oa | Tg) + EG) dr see (A3) 


where p,,= critical density of the substance, and g,(r) is the radial distribution 
function for a liquid at fairly low temperatures. 

This equation is derived from a rather more general model than the simple 
Dieterici liquid, and as it stands it does not lend itself readily either to numerical 
computation or to use in the main body of the present paper. By expressing the 
quantities in it in terms of the a and b of eqn. (A1) we shall derive a formula 
which, while of more restricted applicability, is more readily useful for those 
liquids for which (A 1) is adequate. 

The constant 6 is four times the effective volume of the molecules, whence 


b=4N(47/3)(r9/2)> or ryp=(3b/207N)B (A4) 


where 79 is the effective diameter of a molecule. 
It can easily be shown (see, for example, Slater 1939) that 


a=1N? | Anr2E(r) dr 
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where (r) is the interaction potential energy. Following Lennard-Jones and 
Devonshire (1937) we take E(r) of the form 

E(r) = —p/r>+A/r”. 
(We take the limit ”— oo, instead of n = 12, i.e. we assume perfectly rigid molecules 


of finite size. Fig. 2 shows that the two assumptions differ little.) With this 
expression for E(r) we deduce 


GeiruN yor) ee SRS (A5) 
and thus u is given by 
p= srealéaN*=9abi4r° Ney 4 eas (A6) 
Fowler introduces a quantity «’, defined by 


in PE(r)go(r) dr =a! he PE(r)dr, wasn (7) 


He shows that for a solid close-packed lattice «’=1-:91. In any case «’ >1 and 
as the final expression for P will be proportional to (3«’ + 1)1/4, the extreme values 


Table 2+ 
Substance Per M* Ver b Per a ails eae wg theor, -surf.tens- 
Argon A O53" 39°95 9 7522) 3776 48-6 2:03 1:20 8:29 54:5 52-5 


Nitrogen N, O-31t5 828-0) 90:25 845-4 aoe) 2°05 lori Dow Ose»  _ (LG: 
Methane CH, 0-162>, 16:0 98-8 49:4 46:4 3°32 1-36 WEE) “Uses (BPP) 
@aumnonia NH, 0-235 17:0 72-3 36:1 113-0 4239 1:45 Soll Oo ommmOU [7 


Hydrogen 

chloride HCl 0-42 36:5 86:9 43-4 82:7 4-64 1-46 9-02 72°6 67°8 
Water H,O 0-329 18-0 54:7 27-4 221 49-8 2:66 6:90 108 52-6 
Ethyl ether 

(C,H;),0 0-263 74:0 271 135) 36:0 19°5 Be ies, ANS 212 
Carbon 

tetrachloride 

CCl, 0-558 154 303 £151 45-6 31-1 DEK SEY! WSS) 220 


Benzene C,H, 0:304 78:0 257 128 48:3 23:6 2m Ono's 206. 
Chlorobenzene ; 
C,H,Cl 0°365 112:5 308 154 45:2 30°8 2°36 18:9 261 245. 


+ Units are as follows: per, g/em?; Vy, cm3/mole; b, cm3/mole; per, dyn/cm? x 108; 
a, dyncm*/mole? x 10"; a4, dyn1/* cm/mole?/? x 10°. 


#’=1 and «’=1-91 give values for P differing by only 16%. In fact Fowler’s 
assumption of close-packing is probably not far from the truth, and therefore 
#’ is taken to be 1:91. Then eqn. (A3) reduces to 


aN? M* (3x! +1) 


pts 7 ll 7 cis ee (A8) 

For a cubic close-packed lattice d=(b/2\/2N)"= nen es (A9) 

and for a Dieterici liquid py=M*/V.=M*/2b.  — .wseeee (A 10) 
Substituting (A 9) and (A 10) in (A 8) we obtain 

TIA a ed a es etre (A11) 


where w =[3(3«’ + 1)/647N!]4, ‘Thus, provided we take a suitable value of 
a’ (x =1-9, see above), P is now expressed wholly in terms of known quantities. 
aand bare found from the critical constants, by means of the well-known relations. 
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for a Dieterici fluid: b=4V..=4M*/o,, (where p,, is in g/em® and 6 in 
cm3/mole!), a=4e2b2p,,=29-55b’p,, (where p,. is in dyn/cm? and a in 
dyncm!/mole?) and w =5-87 x 10-3. Table 2 shows values of the various relevant 
quantities for a number of substances and we see that, with the exception of 
water, the values of P deduced from (A 11) are in moderate agreement with those 
deduced from the surface tension by: the definition (A 2). 
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Classical Theory of Compressibility 
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ABSTRACT. By extending Green’s treatment of the classical virial theorem, a formula 
for the bulk modulus of an assembly of atoms, when classical physics applies, is derived. 
‘This formula involves the binary and ternary distribution functions. 


on the theory of liquids, a simple derivation is given of the virial theorem 

for an assembly of atoms to whose interaction classical mechanics applies. 
In this paper the same technique is used to derive a formula, evidently new, for 
the compressibility (or, rather, for the bulk modulus —V(@p/@V),). To 
establish the notation and method to be used we begin by recapitulating Green’s 
treatment of the pressure. ‘The coordinates of the N atoms are combined in a 
single 3.N-vector x, and their interaction potential denoted by ®(x). Let them 
be confined in an enclosure (which we shall vary without changing its shape) 
of volume V. Consider now the integral 


I(f)= J OSG xy a uk: ca ee (1) 


* Now at Department of Chemistry, Duke University, Durham, North Carolina. 


T the appendix to the second paper (Green 1947) of Born and Green’s series 
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taken over the enclosure N times. To ‘differentiate with respect to V’, we replace 
x by Ax and carry out [d/dA....],_,: this is the operation 3V d/dV. Then we have 


3VdI( f)/dV =|Ad/aa ee | 2°™xfAx) ‘abe =3NI(f)+1(x.Vf). oil. (2) 
The Helmholtz potential A is given by the phase integral 
Si Osa(eT Clee pe ea ies (3) 


where a is a constant of no interest here. Differentiating (3) with respect to volume, 
by (2): 
—3V(OA/OV)_ 4" = a(kT)* (3 NRT I(e~ 2!**) — I(x. VOe— 82), oes, (4) 


Now the thermal average of a function f(x) is 


(f= fle @PPY He“ WP2), cane, (5) 
Then, from the definition of the pressure p=—(0A/0V),, eqn. (4) reduces to 
the classical virial theorem 


SPV = 3 NRTA (KLIP) A i FO gees os (6) 
Eqn. (6) takes a more usual form if we assume that 
Dr Ol) eae = at pegs tes eee (7) 


where the r, are the particle coordinates (projections of x), r,,;=",—-¥r,, and 
X,=322,_,; is the sum over pairs of (7,7). If (7) holds, then 


(x . V)"® =24¢,,(7;), where ¢,(7)=(rd/dr)"(r). we eee (8) 
It is then convenient to introduce the multiple probability densities in coordinate 
space, ”,(r,, r,,...",), which we normalize so that J(my)=1 and n,_,= are 
where i ‘stands for | " Pr and for convenience we suppose my made 


symmetrical in the particle coordinates. Then, using (8), eqn. (6) reduces to 


(1,2) 
3pV =3NkT—4N(N-1) | ial Wists) (113) sual eee (9) 
Passing to the derivation of the compressibility formula, we write (6) in the 
form (3pV —3NRT)I(e~ 2**) = — I(x. VO e— 2/**) 


and apply 3Vo/dV again, obtaining 
3V[A(3pV)/OV ]pl(e~ *!"*) + (3pV —3NRT){3NI(e~ 9") — (RT) U(x. Ve #*7)} 
= — {3NI(x.VOe~ #**) + I[{(x. V)?O — (x. VO)?/kT}e~ P77}. 
That is, by using (6), 


9V[O(pV)/OV] p+ (x. V)?®) =(RT)1[ (x. VO)?) — (x. V®)7] wee (10) 
=(k )-"D; say. 
The implications of (10) may be seen by applying (7) and (8). ‘Then 
((x. V2) =1M(N— fey eae ONd.(12) aaa nee (11) 


with the obvious abbreviations. To reduce the ‘fluctuation’ term D of (10) 
we have, by (8), 


((x.V@)?) =$M(N=1) { [°° a(1,2)$,2(12) 
+2(N—2) ["”ng(1, 2, 3)px(12)44(23) 
+3(N-2)(N—3)[[ one (ae 2)4(12)]"} sntoaacanp be eree (12) 
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Combining the square of the final term of (9) with (10), (11) and (12), we obtain, 
finally, . 


9VIA(PV)/OVIn=EN(N—1) { — fo” a(t, 2)64(12) 
+(ATY>Y f°” ng(1, 2)p4°(12) 
+2(N—2) ["" ng(1, 2, 3)4x(12) $4(23) 


—(2N-3)( f°" n9(1, 2):(12)) [p= sae oe (13) 


Thus, subject to approximation (7), the compressibility depends on m, m, and 7s, 
but not on the higher 7,; whereas the pressure depends on m, and 7, only. By 
dropping the essentially positive coefficient of (RT)! in (13) we obtain, for a 
fluid, the inequality 


9[4(PV OV Ip + 4v* | ~ g(r)be(r)d*r>0. 


Here we have set m,(1, 2)=m,(1)m,(2)[g(742) + O(1/N)], v= N/V, and assumed 
N very large. On the other hand, (13) itself simplifies, if we substitute the 
Kirkwood approximation 


(1, 2, 3) =mg(1, 2)n9(2, 3)m9(3, 1)/m,(1)m,(2)m,(3) 
to A(pV)/OV Ir +4? |” g(r) ho(r) dr = (AT) { 40 |” g(r) 2%) er 
+2 [f° e¢rrdg(ra)lg(ra2)—Vdr(ra)pa(ra)d?ry re bee (14) 


If the radial distribution g(r) and the interaction ¢(r) are sufficiently well known 
in a particular case, (14) could provide a test of the accuracy of the Kirkwood 
approximation. When used to derive expressions for the virial coefficients, 
eqn. (14) would also provide a test of Green’s (1947) integral equation for g(r) 
independent of that given by the virial eqn. (9), which was recently put to this 
use by Rushbrooke and Scoins (1951). More general expressions, analogous 
to (13), could clearly be derived for the elastic constants, in the case of a crystal. 
Returning to the general and (for classical physics) exact result (10), the 
‘fluctuation’ expression D clearly must vanish at least as fast as T, when T->0, 
for the formula to be applicable. For a classical crystal D->0 as the temperature 
goes to absolute zero, because the atoms are finally in fixed positions, at low 
temperatures the r.m.s. deviations of the normal coordinates, and hence of the 
distances between atoms, vary as 77”. Hence D may be expected to vary as T, 
and the right-hand side of (10) to tend to a finite value (calculable in terms of the 
normal lattice vibrations) at absolute zero. In reality, the presence of the term 
in D shows that the formal identity of the classical and quantal formulae for 
thermodynamic pressure (Price 1950) cannot hold for compressibility as, owing 
to the zero-point motion, D/kT must tend to infinity at absolute zero. It is hoped 


to publish later the quantum treatment corresponding to the classical theory of 
the present paper. 


(1, 2) 
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ABSTRACT. By measuring the energy of photo-protons from deuterium the x-ray 
spectrum of a betatron, of maximum energy 20 Mev, has been deduced. A small 
difference between the experimental determination and the theoretical curve in the 
10 Mev region is found. 


(1948), Koch and Carter (1950), Wang and Wiener (1950)) have been 

made to determine the x-ray intensity spectrum of betatrons. The simplest 
and most elegant method seems to be due to Wang and Wiener, who used a plate 
soaked in heavy water for a 10 Mev machine. For higher x-ray energies there 
is the added complication due to photo-protons produced in the emulsion from 
(y, Pp) reactions in silver, oxygen and bromine. This difficulty has been resolved 
by using heavy wax film targets. ‘The heavy wax has a chemical formula 
C,,D,.. (n=30) and contains about 25°% of deuterium by weight. ‘The film, 
5 microns thick, was supported on a thin backing of ‘Cellophane’ which was in 
turn fastened to a substantial metal frame, the latter being made so that the 
wax fitted closely to the surface of a quarter-plate C2 Ilford emulsion 200 microns 
thick. 

The emulsion and film were exposed perpendicularly to the axis of the 
collimated x-ray beam at a distance of about 150 cm from the 0-010-in. thick 
tungsten target for a period of two minutes, the x-ray output from the machine 
being about 10 réntgens per minute at one metre, as measured on a thick-walled 
ionization chamber, with a maximum energy of 18-5 Mev. In order to estimate 
any possible background effect an ordinary paraffin wax film was used as a 
target and exposed under similar conditions. ‘The plates were developed in a 
solution of Azol and 1°, potassium bromide. 

Several square centimetres of emulsion were scanned and only photo-protons 
with angles of dip between 5° and 45° were accepted. Over 1,000 photo-protons 
were measured in the heavy wax plates, whilst in the control plates only about 
40 tracks entering from outside the emulsion were found. ‘The majority of the 
latter protons had energies less than 2mMev. The number of photo-protons 
produced by x-rays at different energies from 4 to 18-5 Mev is shown in the 
histogram (fig. 1). The smooth curve is the theoretical distribution of the 
number of protons expected from the Bethe~Peierls theory of the photo- 
disintegration of the deuteron and from the Bethe-Heitler theory of 
bremsstrahlung. The intensity spectrum is given in fig. 2, and it can be seen 
to differ slightly from the theoretical curve. It exhibits similar characteristics 
to the spectrum found by Koch and Carter using pair production, in a cloud 
chamber, to measure the number of quanta. For instance, the increase in the 


S18 previous attempts (Lasich and Riddiford (1947), Bosley et al. 
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number of quanta in the 10 Mev region and the rapid fall towards the high energy 
tip are illustrated. These two factors are obviously dependent upon the point 
of normalization of the theoretical curve and the experimental histogram. If, 
however, the theoretical proton—energy curve is transformed, by adjusting the 
width of the energy intervals into a rectangular distribution, the conclusions will 
be independent of normalization. ‘This transformation has been carried out 
in fig. 3, and the experimental results again illustrate an increase in the medium 
energy regions and a more rapid fall off towards 18 Mev. The causes of these 
two effects may be inherent in the betatron since it is unlikely that the theory of 
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bremsstrahlung is at fault. ‘The increase may be tentatively explained by 
assuming that the electrons, after acceleration, strike the tungsten target and are 
scattered back into the doughnut. After further revolutions the scattered 
electrons hit the target again. The process may occur more than once, and the 
measured spectrum would be the summation of more than one bremsstrahlung 
curve, one due to the primary beam and the other due to lower energy scattered 
electrons. It would be of considerable interest to carry out a spectral analysis 
of some non-cyclic accelerator, as this would rule out the possibility of the 
electrons striking the target more than once. 
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ABSTRACT. Some measurements on the angular distribution of thin target radiation 
produced in a synchrotron are reported using ionization chambers to measure the x-ray 
intensity. Quantitative agreement has been obtained with the bremsstrahlung—multiple 
scattering theories of Schiff and Lawson for the case of a platinum target in the range of 
electron energies of 10-14 Mev. 


§1. INTRODUCTION 

HE theoretical angular distribution of thin target radiation was first formu- 

| lated by Schiff (1946), who attributed the angular spread to multiple elastic 
scattering of the electrons and their subsequent radiative collisions with the 

target nuclei. More recently, an essentially equivalent theory has been presented 
by Lawson (1950), in which the simplification lies in the form taken for the polar 
diagram of a single radiative transition. "The common features of the two theories 
are: (i) the prediction of symmetry about the line of maximum intensity, which 
lies in the direction of the incident electron beam, and (ii) the inverse proportion- 
ality of the angular width and the incident electron energy. In practice, Lawson’s 
formula is more directly useful since it avoids the divergence of Schiff’s formula 
for small values of the angle 6. Because of this divergence, Schiff’s formula is 
applicable only to angles somewhat greater than ,./E, where p is the rest energy of 
the electron, and £ the total electronenergy. We have found it possible to remove 
this limitation by a simple modification of the Schiff formula which makes it 
applicable for all angles. This was achieved by using the asymptotic property of 
the exponential integral, namely Ei( —,)—Ei(—‘y2)—>I1n(y4/ye) as ¥4, v2->0. 
The second exponential integral term is determined from the condition /(@)/J(0)—1 
as@—0. This modification then allows for the finite width of the radiation polar 
diagram, which plays an important role at small angles. ‘The expression thus 
obtained is analogous to that of Lawson. Calculations made for a tungsten 
target, using the modified formula, differ by less than one per cent from the 
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published curves of Schiff. The following are the expressions of the radiative 
polar diagrams for the three cases: 

16) —Ei{—(£6)" In (1832Z-19)/1510-82} Schiff 1 
1(0) In {1510-8t/p2 In (183Z25)} — 0-5772 (Sebi) eae ae 
1(@)  —Eif{—(E6)? In (183Z-1)/1510-8¢} + Ei{ — (£8)?/1-787 7} 


1(0) In {1510-8¢/p7 In (183Z-"8)} — 0-5772 

(Modified Schulyy ur . oe (2) 
1(6)  —Eif—(£6)?/(440¢ + 2u?)} + Eif — (£6)?/2u.7} 4 
70) ~ Sacer (E717 iolssirs «a (Lawson)j 0°70 (3) 


where J(6) is the intensity, for a given target distance, at an angle 6 with respect to 
the direction of incidence of the electrons, E the electron energy in Mev, » the 
rest energy of the electron in the same units, and ¢ the target thickness expressed in 
radiation lengths (Rossi and Griesen 1941). These formulae are applicable to: 
target thicknesses in the range from 10-* to 10~! radiation lengths. 

Comparison with theory of the experimental results obtained with a betatron 
or synchrotron is complicated by a lack of knowledge of the exact conditions under 
which the electrons strike the target. Radial oscillations of the electrons contri- 
bute the main uncertainty in the knowledge of their angle of incidence, so that the: 
effective thickness of the target may be slightly different from its geometrical 
thickness. However, this factor is usually of secondary importance since the 
intensity at an angle @ is relatively insensitive to small changes in the target: 
thickness. 

Approximate agreement with theory has been reported in several experiments 
(Lawson 1950, Lees and Metcalfe 1950, Allen-Williams and Appleyard 1949 a, b) 
but in most cases, the target conditions have made a quantitative comparison 
impossible. Preliminary reports (Baldwin, Boley and Pollock 1950, Lanzl and! 
Hanson 1950) confirm, with reservations, Schiff’s theory, although deviations from 
theory in the case of heavy elements were found by Baldwin et al. It is the 
purpose of this paper to present data on the angular distribution obtained with a 
14 Mev synchrotron using air-ionization chambers for the intensity measurements.. 


§2. EXPERIMENTAL ARRANGEMENT 
(i) Production of the X-ray Beam 

The accelerated electrons were made to spiral slowly out from the equilibrium: 
orbit at the end of the acceleration time and the x-ray beam was produced in a 
platinum target, 0-005 in. (0-045 radiation length) thick, which was mounted out 
from the back of the injector anode perpendicular to the electron orbit. The target 
was a flag three millimetres square, which ensured that the x-ray beam was pro- 
duced well clear of the injector anode. ‘The azimuthal position of the target was 
chosen in order to reduce scattering from the magnet poles and energizing coils 
toaminimum. ‘The peak electron energy was obtained from the experimentally 
determined magnetic field characteristics. Using this magnetic calibration the 
gamma-neutron threshold in silver has been determined and the value obtained 
is in good agreement with the value given by Baldwin and Koch (1945). 


(11) Detection and Monitoring 
The x-ray intensity was measured by a cylindrical aluminium ionization. 
chamber of one centimetre wall thickness. This instrument, which is direct 
reading, was shown to be linear for the intensities encountered during the experi- 
ment. Good resolution was obtained by using the chamber in the end-on position 
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at a distance of 3} metres from the target; under these conditions, the chamber 
subtended an angle at the target of about one degree. A thin-walled integrating 
chamber, used for monitoring the x-ray output was placed in a standard position 
and in such a way as not to interfere with the beam proper. To allow for fluctu- 
ations, the intensity at each angular position of the ionization chamber was read 
regularly during a suitable period and the integrated intensity estimated. The 
latter was normalized to a standard irradiation using the corresponding monitor 
readings. For some runs, the two chambers were interchanged and no significant 
difference was found for the results of these runs. 


(iii) Scattered Radiation - 


The effect of scattered radiation was estimated by placing a lead absorber, 
seven inches thick, directly in front of the ionization chamber, which was mounted 
ona light frame to minimize the effect of scattering. ‘The measured intensity under 
these conditions was entirely due to scattered radiation and was comparable with 
the scattered intensity present in the absence of the lead absorber. The residual 
intensity thus measured was found to be distributed isotropically in the region of 
the beam, and never exceeded one-half of one per cent of the x-ray intensity at the 
beam centre. ‘The contribution of scattered radiation is therefore negligible for 
the present experiments. 


§3. RESULTS 


Using the inverse square law, a small correction was applied in order to reduce 
all readings to the same target-chamber distance. Each set of readings was given 
a preliminary examination by plotting the intensity as a function of the arbitrary 
angular position of the ionization chamber. In every case the distribution was 
found to be completely symmetrical. ‘This procedure is necessary since the centre 


Table 1. Characteristic Angular Width « as a Function of the Electron Energy EF 


Eu Plane of angular 
BABY PAGES) rey. deg.) distribution 

10-2 VBS) 80:6 Horizontal 
12-3 6-2 76°3 os 

12-7 5:9 75-0 = 

13-1 6-0 78:6 vs 

13°8 5°5 75°9 a 

11-6 6:7 Tika] Vertical 
12-1 6:33 75°8 nf 

12-3 6-1 81:8 

Mean 77:8 


of the beam was observed to alter in the case of distributions measured in the hori- 
zontal plane through the target as the electron energy was altered, as predicted by 
Allen-Williams and Appleyard (1949a). From these curves the characteristic 
angular width « subtended at the target between the lines of maximum and half- 
maximum intensity are obtained; the results are tabulated in table 1. It can be 
seen that the product of angular width by the peak electron energy is, within 
experimental error, a constant characteristic of the target, as predicted. ‘The value 
of this product as given by Schiff’s formula for a target of thickness 0-045 radiation 
length is 76-7, as compared with our experimental mean value of 77-8. _ Lawson’s 
theory gives 80-9. From this it follows that there is no significant difference 
between the effective and the geometrical thickness of the target. 
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An overall test of the theory is possible by compounding the various runs first 
by converting the independent variable 0, measured with respect to the beam centre, 
to the product of angle and peak energy E@. The results are given in table 2 and 
are displayed in the figure for ease of comparison with theory. ‘The errors shown 
are the probable errors obtained from table 2 assuming the spread of values to be 
entirely due to experimental errors. Comparison is made with (i) the Schiff 
curve, (ii) the Lawson curve, and (iii) the modified Schiff curve, each being 


Table 2. Experimental Values of 100J(@)/J(0) as a Function of E@ for a 
Series of Electron Energies F (Mev) 


E 10:2° “42-3 «1297 1331 13s 6 412 13 ee 
E6 Hd HOY HG) Hw HO VO Voew eo 
10 97:0 98-5 96:3 94-7 96:0 98-0 97:0 97-0 96:8+0:3 
20 90-7 95:8 90:3 87:0 87:38 94:0 91:7 94-2 91-4+0-7 
30 82-3. 906 83-0 79:7 80-7 88-0 85-0 -89-5 84-8+1-0 
40 740 83:2. -75-0-; 72:7. 373-65, 185 “777 eso 77-3411 
50 67:0 75-3 67:8. -65-6 66:5  70°8- 69-6 76-0 69:8+0-9 
60 60-8 65:1 60-2 59-6 59:3 623 62-0 68-0 62:2+0-7 
70 55°55. 56:0 53-1 54:3 53:0 * 55-0 54-3 60-0 55-1+0°5 
80 50-4 49-1 47-3. 49-1 48-0 - 48-2 47-6 «51-5 48-9 + 0-4 
90 43-8 44:0 42:7 44:3 43-5 41-3 41-8 44-0 43-2403 


H (1)—H (5) distributions in horizontal plane through target. 
V (1)-V (3) distributions in vertical plane through target. 


----— Schiff 
Lawson 
Modified Schiff 


Experimental Points 


1002(@)/7(0) 


EO (Mev. deg) 
Radiation polar diagrams for target thickness of 0-045 radiation length. 


calculated under the conditions mentioned above. Although the precision of the 
experiments is not sufficient for a decision to be made between the forms of the 
distribution, there is slightly better overall agreement with the modified Schiff 
curve than with the other two. The deviation at small angles for the Schiff law is 
due to the breakdown of the theoretical assumptions for angles less than 1/2E 
radians. 

In some sets of measurements, two peaks were observed in the angular distri- 
bution curves. ‘This phenomenon, which has been found for several different 
electron energies, was due to two separate beams as shown by the existence of two 
sets of Geiger—Miiller counter pulses separated by about one millisecond. There 
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were found to be two different acceptance times for electron injection, one of which 
gave rise to two beams and the other to a single beam. Consequently care was 
exercised in choosing the appropriate injection conditions during the course of the 
experiment. Double beam distributions were excluded from the data presented 
above. 

§4. CONCLUSION 

As noted above, all distributions were found to be completely symmetrical, in 
contrast with the results of Lees and Metcalfe (1950) and Allen-Williams and 
Appleyard (1949 b). In our case asymmetrical distributions were not observed 
after precautions were taken to eliminate scattering from dense bodies near the 
synchrotron. Also the use of an external target removes more completely the 
effect of scattering from the magnet poles and the small yet complicating effect of 
the residual betatron beam. 

The formulae of Lawson and Schiff (modified form), as applied to a target of 
platinum of 0-045 radiation length, are in agreement, within the statistical error, 
with the experimental results given above. Brief details of a similar experiment 
by Lanzl and Hanson (1950) have recently appeared, in which angles of half 
intensity for 17 Mev electrons impinging on gold targets of various thicknesses were 
slightly less than predicted by Schiff’s theory. On the other hand, Baldwin, 
Boley and Pollock (1950) found good agreement for targets of low atomic number 
but deviations from theory, which increased with atomic number of the target 
material, for thicknesses of the order of 10-? radiation length. ‘The only known 
difference in experimental conditions was the use of electrons of energy 70 Mev, 
and therefore a possible explanation of the deviation found in this case lies in the 
increasing importance of the bremsstrahlung process in determining the angular 
distribution at an energy as high as 70mMev. From the figure it is possible to 
suggest that even in the 10-14 mev range, the bremsstrahlung distribution, which 
plays a major part at small values of 0, is somewhat narrower than assumed by 
either Schiff or Lawson. For lower or higher energies, deviations may be 
expected, in the former case on theoretical grounds (Schiff 1946), and in the latter 
in the light of the experimental work discussed above. 

Further work seems to be called for in order to decide at what energies the 
assumptions of the Schiff and Lawson theories regarding the bremsstrahlung polar 
diagram become inadequate. 
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Disintegration Protons from the Bombardment of Beryllium 
by 7°7 mev Deuterons 
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ABSTRACT. The deuteron beam from the Liverpool 37 in. cyclotron was used to 
bombard a thin beryllium foil in a vacuum chamber. The photographic plate method 
was employed in detecting the charged particles emitted from the reaction. The angular 
distribution of the elastically scattered deuterons shows a secondary maximum at 65°. 
Angular distributions of the two longest range proton groups are compared with the 
recent theories of the angular distributions of the (d, p) and (d, n) reactions. ‘The results 
indicate that both the ground state and first excited state of the residual nucleus !°Be have 
€ven parity and spin either 0, 1, 2 or 3. 


§1. INTRODUCTION 

TTENTION has been recently given to the study of the distribution of 
A particles emitted in nuclear reactions with respect to their angle of 
emission. ‘The investigation of such distributions for protons resulting 
from the bombardment by deuterons of about 8 Mev indicates interesting 
characteristics at small angles of emission. Pronounced maxima are observed in 
the cases of oxygen (Burrows, Gibson and Rotblat 1950) and aluminium (Holt 
and Young 1950, Gove 1951) and are attributed to the possibility of the stripping 
of the deuteron in its passage past the target nucleus. In the light of the recent 
theories proposed by Butler (1950) and Bhatia and Huang (unpublished) such 
observation can give information regarding the spin and parity of the 

appropriate state of the final nucleus. 

In the present work the angular distributions of the two longest range 
proton groups from the (d,p) reaction with beryllium are investigated. Also 
the angular distribution of the elastically scattered deuterons from this mono- 
isotopic light element is found to show features similar to those observed by 
Guggenheimer, Heitler and Powell (1947) in the scattering of 6-5 Mev deuterons 
by oxygen and nitrogen and in the scattering of 8 Mev deuterons by neon 
(Middleton. and Tai 1951). 

§2. METHOD 

A thin beryllium foil 0-325 Mg/cm? thick was bombarded in a vacuum 
chamber by a beam of deuterons of 7-70 Mev mean energy provided by the 
Liverpool 37 in. cyclotron. ‘The cross section of the beam had a diameter of 
1:5 mm and its angular divergence was about one-third of a degree. The 
emitted particles were detected in Ilford C2 emulsion 200 u thick and examined 
with a binocular microscope with magnification of about 900. <A detailed 
description of the apparatus used in the present experiment, as well as the method 
of measurement, was reported by Rotblat, Burrows and Powell (1951) and 
El-Bedewi (1951). 

The experiment was carried out with two different exposures to allow for the 
expected sharp rise of the emitted particles near the forward direction. 


Deuteron Bombardment of Beryllium ) 65 


A normalization was obtained by measuring the intensities of different groups of 
particles at an angle of 30° in both exposures. The mean ratio of the intensities 
of the forward to the backward exposures was found to be 0:51+0-01. Also, 
for the investigation of particles emitted within a continuous angular range the 
foil was tilted to the direction of the incident deuterons by either 45° or 135° to 
suit the forward or backward exposures respectively. 

In examining the plates, a certain area of the emulsion surface was carefully 
explored and all the tracks belonging to the elastically scattered deuterons and 
the two longest range proton groups from the reaction °Be(d, p)!°Be were 
recorded. ‘The position of measurement was chosen such that the angles of 
approach of the particles entering into the emulsion were between 6° and 10°. 
The loss of tracks through scattering out of the top surface of the emulsion was 
then found to be less than 2% everywhere. ‘Throughout the whole series of 
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Fig. 1. Histogram of disintegration products Fig. 2. Angular distribution of proton group 1 
of beryllium bombarded by 7:7 Mev in the centre-of-mass system. 
deuterons at 70°. 


measurements only three tracks were observed to run into the glass backing of 
the emulsion, as might be expected from a comparison of the ranges of tracks 
under consideration with the thickness of the emulsion used. Since the 
percentage of tracks scattered out of the emulsion was small, the variation of this 
loss at various angles of observation should be even smaller. ‘Therefore, in the 
results presented below, no correction of this loss was made. ‘The errors shown 
in the figs. 2, 3 and 4 are those of a statistical nature only. 


§3. ANGULAR DISTRIBUTION OF THE PROTON GROUPS 
A representative histogram for the upper end of the particle spectrum at an 
angle of observation of 70° is shown in fig. 1. Groups 1 and 4 are identified, 
using known Q-values, as being due to protons leading to the ground state and 
to a level at 3-375 mev of the residual nucleus }°Be. Groups 2 and 3 are due to 
oxygen contamination in the target. The angular distributions of groups 1 and 4 
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have been measured by a method similar to that described by Chadwick, May, 
Pickavance and Powell (1944). The number of tracks belonging to group 1 is 
known directly from the plotted histograms since it is well separated from all 
other groups. In the case of group 4 great care has been taken to discriminate 
that group from the neighbouring tracks resulting from the foil contamination. 
By fitting Gaussian curves, reliable counting was achieved at all angles of observa- 
tion except for those between 30° and 60°, where there was an overlapping with 
group 3. Fortunately, the angular distributions for protons from the reaction 
16Q(d, p)!7O have been worked out in this laboratory (Burrows et al. 1950) under 
the same conditions, and one can calculate the ratio of the numbers of protons 
belonging to the ground state to those belonging to the first excited state of 
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of-mass system. 


Q at a certain angle of emission. Since the ground state group 2 was well 
separated in the spectrum from the beryllium target, it was possible, by making 
use of this known ratio, to estimate the number of protons due to contamination 
in group 4. This procedure was repeated for each angle where the two groups. 
overlapped. 

After calculating the relative intensities of each group at the various angles 
of observation, they were normalized, transformed to the centre-of-mass system 
of coordinates and then represented as shown in figs. 2(a) and 3(a). These 
results are based on the measurements of a total of 1664 tracks for the 
group 1 and over 22 different angles and 3846 tracks for the group 4 over 
23 angles. ‘The distribution for the former group shows a forward maximum 
at 25° with a smaller one at 75° and a slight rise at the backward angles. The 
other group appears to show a forward maximum at about 20°, but the statistics 
are not good enough to establish this definitely. 
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§4. ANGULAR DISTRIBUTION OF THE ELASTICALLY SCATTERED 
DEUTERON GROUP 

In order to count the number of particles belonging to the elastically scattered 
deuterons from beryllium only, histograms for the accepted tracks in the 
neighbourhood of its range are plotted for each angle of measurement and then 
fitted with Gaussian curves. The angular distribution was then calculated as 
previously described and represented as shown in fig. 4, where the scale of the 
coordinates is chosen to be the same as for the proton groups. ‘These results 
are based on the measurements of a total of 6075 tracks over 16 different angles 
at intervals of 25° at the secondary maximum appearing about 60°, and 10° for 
the rest of the distribution. 


§5. DISCUSSION 


The angular distributions of the proton groups show a pronounced increase 
in intensity at some small angles or in the forward direction of emission similar 
to those observed in other reactions (Holt and Young 1950, Burrows et al. 1950, 
Gove 1951). This small angle structure can be interpreted by the mechanism 
recently proposed by Butler (1950) and Bhatia and Huang (unpublished), in 
which they assume that in this region the reaction proceeds mainly through a 
stripping process. Applying the latter theory, Huby and Newns have calculated 
theoretical calculations of the (d,p) reaction with beryllium. The radius of 
interaction used in the theoretical formula was 5-41 x 10- cm. In comparing 
the theory with the experiment the best agreement was achieved only when the 
orbital angular momentum /, brought to the nucleus by the absorbed neutron 
was equal to unity in the case of both of these proton groups. Any other value 
of /, would fail to agree with the experiment unless an impossible value for the 
radius of interaction was used. The theoretical curves, together with the 
experimental points for the proton groups leaving the residual nucleus in the 
ground and first excited states, are shown in figs. 2 (b) and 3 (6) respectively. 
The theoretical curves were normalized so as to fit the experimental results in 
the neighbourhood of 25° to 40°, where the experimental distributions show large 
variations of intensities with respect to angles of emission. ‘This was so chosen 
because the theoretical distributions for different values of /, show marked 
differences in the steepest slopes along the distribution curves. Similar results 
can also be obtained by applying Butler’s theory and assuming the radius of 
interaction to be 4:8 x10-!%cm. In all the measured angular distributions the 
points corresponding to the two smallest angles (174° and 20°) were found to be 
lower than would be expected from the theoretical curves. ‘This may be due to 
loss of particles in the apparatus, since these angles are near to the geometrical 
cut-off of the apparatus, which was known to be at approximately 16°. 

Information about the spins and parities of the appropriate states of the 
final nucleus can be thus obtained if these properties are known for the ground 
state of the initial nucleus. Recent work reported by Schuster and Pake (1951) 
gives the spin of the ground state of *Be to be 3/2. ‘The parity of this state is known 
to be odd. The fact that /,=1 for both the ground and the first excited states of 
10Be indicates that these have even parity and spin either 0, 1, 2 or 3. Hornyak 
et al. (1950) give spin zero and even parity for the ground state of Be. Also a 
gamma transition from that excited state to the ground state is reported by 


Rasmussen, Hornyak and Lauritsen (1949). If the assignment of spin zero, 
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even parity to the ground state is correct, this would exclude spin zero for the 
excited state, but the other possibilities for the spin remain. ; 

With regard to the secondary maximum appearing in the distribution of the 
elastically scattered deuterons, a theoretical estimate of the angle at which this 
secondary peak ought to appear has been worked out by Huby (private 
communication), on the assumption that at this angle the coulomb scattering 
is negligible compared with the nuclear scattering. Using the Born approximation 
with a spherically symmetric nuclear potential effective at a radius of 
5-41 x 10-18 cm, the angle found is 64°, which is in good agreement with the 
observed angle. 
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Ionization in Oxygen by »-Mesons 


We report here preliminary results of measurements of the ionization by jx-mesons in 
oxygen as a function of meson momentum. These show the expected logarithmic increase 
of ionization in the relativistic domain, in agreement with the recent measurements of 
Becker et al. (in process of publication) in krypton, and further give some indication of 
the onset of the polarization effects treated theoretically by Fermi (1940), Wick (1943) and 
Halpern and Hall (1948), which have already been observed in condensed material by Bowen 
and Roser (1951). 

The region of high momentum (p>2 x 101° ev/c), where the polarization effect in a 
gas at normal pressure may be expected to become large, is the subject of work still in 
progress. 

The present data are based on drop counts in a cloud chamber along the trajectory of 
p-mesons, the momenta of which are determined in the Manchester magnetic spectrograph 
(see Hyams et al. 1950). The detailed precautions necessary in such measurements will 
be discussed elsewhere; briefly, the counts were made in oxygen and a water—alcohol 
condensant, maintained before expansion at atmospheric pressure and at a steady 
temperature of 20+1°c, on the positive column of field separated tracks for which the 
condensation efficiency on the negative ions was high enough to ensure complete 
condensation on the positive ions. The results given in the table are based on the exclusion 
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of large energy transfers leading to clusters of more than 30 ion-pairs, but other criteria 
of exclusion have been considered and the results are not sensitive to the precise procedure 
adopted. The standard error of the specific ionization, given in columns (5) and (6) of 
the table, is that obtained directly from the spread of values among the various tracks of 
each momentum range measurement. This spread is in good agreement with that 
calculated on the assumption that the distribution of secondary ionization corresponds to: 
Rutherford scattering in the primary encounters. 
The group of lowest momentum tabulated falls outside the efficient range of 
' spectrograph; the particles used in this group were therefore selected as having a range 
in lead between 15 and 35 cm. 


Table 
(1) (2) (3) (4) (5) (6) 
3-6 x 108 33 474-5 19,028 40-4-L0-8 44-9-10-9 
6-10 108 9 140-2 6,237 44-44-1-0 49-641-1 
1-2 x 10° 22 329-1 14/493 44-741-2 49-9413 
2-3 x 10° 16 241:°5 Aes 47-6+1:2 522823 
3-4 « 10° 22 356:2 17,774 49-9-+0-9 55°6-+1-0 
4-10 x 10° 14 235°1 12,427 53:0-+1:3 58-8+1-4 
~1-5 x 101° G2 1961 10,327 52°8+1-1 58-6--1:2 
~3x 10° 5 85:1 4,641 55-142°8 60-9+-3-1 


(1) Momentum range (ev/c); (2) number of tracks; (3) total track length (cm); (4) total number 
of drops on positive ions; (5) mean ionization per cm (observed); (6) mean ionization per cm 
(reduced to pure O,, S.T.P.). 


In the table the last column gives the ionization reduced to oxygen at standard pressure 
and temperature and also includes a correction (of the order of +4°/) for the overlapping 
of the drop images. 

The experimental values are compared with those given by the Bloch (1933) relation 
in the figure, where the experimental points and the theoretical curve are fitted at 
210° ev/c. The normalization of the curves yields a mean energy expenditure per 
ion-pair produced of 31:0+0-8 ev.. There is no significant discrepancy between the form 
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of the theoretical curve and the observations below about 10! ev/c; the last two 
experimental points, which refer to meson groups of higher mean momentum, appear, 
however, to fall significantly below the Bloch curve. A crude calculation to the 
approximation used by Fermi shows that the onset of polarization should take place at 


about 101° ev/c. 
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On Vibrational Transition Probabilities in Electronic 
Transitions 


In a recent paper in this journal, Pillow (1951) has shown how distorted harmonic . 
oscillator wave functions may be utilized to obtain vibrational transition probabilities in 
electronic transitions. According to Pillow, the vibrational transition probabilities obtained 
from these distorted harmonic oscillator wave functions are in good agreement, at least — 
for low vibrational levels and for molecules with small anharmonicity, with those calculated 
from the appropriate Morse anharmonic wave function. The advantage claimed for this 
method is that it obviates the large amount of tedious calculation which is required to obtain 
the exact Morse wave function. 

There is, however, another point which should not be overlooked in all these calculations. 
In computing.the vibrational transition probabilities from the overlap integral of the 
vibrational wave functions, one makes the tacit assumption that the electronic moment 
integral R,=J.’/My,"dr is independent of the internuclear distance r (Condon 1928). 
That this is not the case has been shown by explicit calculations of electronic transition 
probabilities and oscillator strengths as a function of internuclear distance (Coolidge, 
James and Present 1936, James and Coolidge 1939, Mulliken and Rieke 1941, Bates 1951). 
Whether this variation of R, with r has any marked influence on the value of the vibrational 
transition probabilities has not yet been investigated in any detail. The only work along 
these lines of which the author is aware (Shuler 1950) seems to indicate that, at least 
for OH? X+—?ID), the variation of R, with r is a more important factor in the calculation 
of vibrational transition probabilities than is the use of more exact wave functions. 

Until some additional information is available on this point, or unless one can take 
explicit account of this variation of R, with x in a particular case, one probably should not 
place too much reliance on the numerical results of any calculation of vibrational transition 
probabilities in electronic transitions, regardless of how well the mechanical anharmonicity 
may have been taken into account. 


Applied Physics Laboratory, Kurt E. SHULER. 
The Johns Hopkins University, 
Silver Spring, Maryland, U.S.A. 
9th October 1951. 
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The Spectrum of NS 


A rotational analysis of the y-system of NS, which was first observed by Fowler and 
Bakker (1932), has recently been given by Zeeman (1951). We had independently carried 
out an analysis of the same bands, though from spectrograms of much lower dispersion 
than used by Zeeman (Hilger large quartz Littrow instrument as against the 6th, 8th and 
9th orders of a 21-ft. grating). Our results are in good agreement with his, and it is perhaps 
not without interest to note that the value of our primary constant, By, namely 0°8283, 
differs from Zeeman’s value, 0°8267), by only Oe. 

However, apart from confirming Zeeman’s results, we can add a little to the published 
data on this molecule. As already reported (Barrow 1950), the 0,0 and 1,0 bands of this 
system have been obtained in absorption. This seems to be the first time that any band 
with v’>0 in the y-system has been recorded. "The measurements are given in the table. 
The values of AG’), are 1386.. and 1387.,cm-!. A very rough value of Xo’We can be 
obtained by assuming that the value of x,’ is the same as for the ground state. This ‘gives 
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Xe'We’ ~ 8.75 and we’ ~ 1403.; cm~!. It may be noted that this new value of w,’ gives, 
with D=4B*/w?, D’=1-16 x 10~-*, in poor agreement with Zeeman’s experimental figure 
of 0:77 10-®. The use of the value D’=1-16 x 10~-® together with Zeeman’s values of 
A;F’(K) would probably lead to a slightly improved determination of Bj, but the change 
would be very small, and it has not been thought worth while to do the calculations. 

The facts about the other known system of NS, the f-system, are, in contrast, rather 
obscure. The transition has generally been considered to be 21I—II by analogy with the 
B-system of NO. However, all our attempts at a rotational analysis from the low-dispersion 
plates on this basis have failed. But it has been possible to produce a partial analysis of 
the 2,0, 1,0, 0,0 and 0,1 bands on the assumption that each sub-band possesses a strong 
Q-branch. Confirmation comes from the agreement of the combination differences, A, F’(J) 
as obtained from these bands with those got from Zeeman’s measurements of the y-system. 


‘Table. Band-Heads in the Absorption Spectrum of NS 


Head Intensity Nair (A) Pyel(GIa-3) 
0,0 &P. 1 2317-27 43139-1 
2 3 15-98 165-0 
P, 2 05-28 365-3 
Q, 3 04-16 386-4 
1 tOu Ps 2 2243-91 44551-2 
Qy 2 OZ is 9 


Note: These are the only bands of NS so far observed in absorption. 


This, together with the fact that the apparent doublet separation in the /-bands differs 
considerably from that in the bands of the ?X+—?IT system, leads us to conclude that the 
transition may be 7A—*II rather than 7JI-?IJ. Preliminary values of B for the upper state 
are : B,=0-694,;—0-005,(v’ + 4). 

Several other features of the ultra-violet spectrum of NS invite inquiry. First, the 
existence of perturbations in the upper state of the f-bands is suggested both by the anomalous 
behaviour of its vibrational constants (for one sub-state, G’ =940-4.u’ —4-8u’?; for the other, 
G’=962-0u’—8:-65u (Zeeman 1951)) and by the appearance of some of the unanalysed 
rotational structure; secondly, there are a number of strong bands observed both by Fowler 
and Bakker and by ourselves which have not yet been assigned. Further work with high 
dispersion should resolve these questions. 


Physical Chemistry Laboratory, R. F 
Oxford University. A. R. DowNiE. 
20th October 1951. Rake 
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On 0 Values in the Resistance of Metals 


Blackman (1951) has recently pointed out that the @ values obtained from electrical 
resistance data are not directly related to the values obtained from specific heat data, since 
only longitudinal phonons interact with the electrons, at least at low temperatures. He has 
calculated values 6;, relating to that part of the specific heat which arises from the 
longitudinal phonons, that is in calculating 61, he used the velocity of longitudinal waves, 
which in turn is obtained from the elastic constants. 

However, an additional complication occurs in that the values of @ derived from the 
elastic constants relate only to the low-frequency phonon spectrum, and are related to the 
specific heat at low temperatures (T* region). The value of 6, associated with the 
maximum frequency in the Debye approximation (i.e. assuming the existence of a sharp 
cut-off frequency) and determining the specific heat at higher temperatures, is somewhat 
lower owing to the decrease of the phase velocity of the phonons with increasing frequency. 
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In the case of a linear chain, the ratio of the high-frequency velocity to the low-frequency 
velocity, and hence the ratio of the two @ values, would be 7/2~1°57, and for a 
three-dimensional lattice the value would probably be similar. 

This dispersion of the phonon phase velocity is verified in the thermal conductivity of 
dielectric solids, treated elsewhere by the author (Klemens 1951). The 6 value of quartz 
derived from the low temperature specific heat is 290° k. The value of @ which has to be 
substituted into the formula for the resistance due to Umklapp processes in order to fit 
the observed temperature dependence at liquid-hydrogen temperatures is about 170° k, 
for in this context 0 refers to the maximum frequency of the transverse phonons. 

In the case of electrical resistance, 8 is determined by the relation between the resistance 
at temperatures T, greater than 6, and T, very much less than 0. At high temperatures 
the resistance is proportional to the total number of phonons present, i.e. 


OIL x2 dx 
@ (es) 


At low temperatures the resistance is due to a diffusion of electrons around the Fermi sphere 
in small steps proportional to w?. Hence 


Re I ae ONE seal coe a (1) 
0 


o ee “Ad: 
Ree | * N(w)ot doce 78 ™ enn Te (2) 
0 (ia lec—ad) 
Both forms can be combined, when taking note of the appropriate multiplicative constant, 
by means of the Bloch—Griineisen interpolation formula 
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If Oy is the value obtained by fitting the resistance at high and low temperatures to (4), 
then, considering (1), we get 6g =(0,,0y)!'2 = 0,,(2/7)!/?=0-86,. It can be seen from the 
values tabulated that this explains, at least in part, the disagreement pointed out by 
Blackman between 67, and OR. 


Metal Na ea Cu 
65, 260 530 505 
Op 202 260 330 
0:8 Oy, 208 425 404 


There is considerable uncertainty in the value of 0),, since the relation 6y—0-8 @, 
applies only to a linear chain. Consequently one cannot expect close agreement. But 
in the case of Li, where the disagreement is considerable, the electrical resistance cannot 
be described by means of a single 8g (MacDonald and Mendelssohn 1950), from which it 
follows that the Debye approximation is not a good description of the phonon spectrum. 
At higher temperatures one must choose 0p ~360° K to fit the observed resistance of Li. 

The above considerations also explain why the @ value derived from the low-temperature 
specific heat, which is always less than 0), agrees often with 0g, which must also be less 
than @;. But it is clear that in view of the many uncertainties involved it will be difficult 
to correlate the @ values obtained from the specific heat and from the resistivity. 


Division of Physics, P. G. KLEMENS. 
Commonwealth Scientific and Industrial Research Organization, 
Sydney, Australia. 
1st November 1951. 
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Neutrons from the Spontaneous Fission of Thorium 


The purpose of this letter is to resolve a point regarding the number of neutrons emitted 
in the spontaneous fission of thorium. 'The early work of Pose (1943) as reported by 
Segré (1951) appears to give a value for v, the average number of neutrons emitted per 
fission, of between 5 and 6. This is considerably greater than that obtained, for example, 
in the slow neutron induced fission of ?°°U for which v is 2-5+0-1 (Atomics 1951). 

The neutron emission from the spontaneous fission of uranium has béen observed by 
several authors (Segré 1951). Littler (1952) obtained a value of 59-5+3-3 neutrons per 
gramme per hour for uranium, and coupling this with Segré’s value of the spontaneous fission 
rate of natural uranium (due to *°8U content) of 24-8+0-9 per gramme per hour we obtain 
vp=2:4+0-2. This agrees well with the value of »y for neutron induced fission. 

It seemed worth while to redetermine the number of neutrons emitted in the spontaneous 
fission of thorium to see if, in fact, the high figure for v of Pose was correct. A slow neutron 
detector had been built similar to, but less sensitive than, one suggested by M. Goldhaber 
(private communication). Our apparatus consisted of a cube (of side 70 cm) of paraffin 
wax with a central cylindrical hole. Eight proportional counters, filled with 1°BF,; and 
coupled in parallel, were embedded symmetrically in the wax around the axis of the hole 
and the output pulses from these counters were fed to a head and main amplifier and then 
into a discriminator and scaler. The absolute sensitivity of the arrangement for detecting 
neutrons from a calibrated radium—beryllium source was about 1°94. However, an accurate 
calibration of the sensitivity was unnecessary since, in this experiment, we simply compare 
the neutron emissions from known masses of uranium and thorium metal under identical 
geometrical conditions. The only assumption made is that the energy spectrum of the 
neutrons emitted in the spontaneous fission of the two nuclides is approximately the same. 
A cadmium cylinder was placed in the central hole of the wax and into this, in turn, were 
put the same number of uranium and thorium rods. The purpose of the cadmium was to 
eliminate any effects which might be caused by slow neutrons re-entering the uranium and 
causing neutron induced fissions. (With thorium such an effect is not possible.) 

The apparatus was housed in a laboratory which was remote from neutron sources and 
the background counting rate of the neutron counters was low, being about 1 count per 
minute for all the eight counters. The counting rates obtained with uranium and thorium 
are given in the following table : 


Metal present Background 
Total Time Time True counts Count/kg of 
Wt. of metal counts* (min.) Counts (min.) per min. metal 
U 33-6kg 4567 ey ML { 303 44:5 9:03 +0-13 
Tage 1s S0S.afe ee. Oe 1-42+0-:09  0-060-£0-004 


* including background. 


From this we see that the ratio of the neutron emissions from uranium and thorium is 
153+10, and since the ratio of the spontaneous fission rates of uranium and thorium is 
(Segré 1951) 164+ 13 we can say that the value of v for thorium is, within the statistics, the 
same as that for uranium, i.e. vpy,/vy=1:07+0-10. This confirms then that Pose’s figure 
of v for thorium is too high. 

We have satisfied ourselves that the contribution of neutrons due to (a, n) reactions with 
light element impurities in the uranium and thorium are negligible in this experiment. 
It is possible that Pose’s determination may have been incorrect because of some impurities 
in the thorium used. For example, a uranium content about 1% by preighie could account 
for the observed neutron emission. 


Atomic Energy Research Establishment, F. R. Barcuay. 
Harwell, Didcot, Berks. W. GALBRAITH. 
14th November 1951. W. J. WHITEHOUSE. 
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Comparison of Geiger and Proportional Counters for 
Intensity Measurements in X-ray Diffraction 


In a tecent investigation with a counter spectrometer’ of the x-ray scattering 
from proteins, the need arose for as pure a monochromatic x-ray beam as possible 
without either the severe loss of available intensity entailed by the use of a plane 
crystal monochromator or the extra labour involved in a double scan with balanced 
filters. In particular, an estimate was required of that fraction of the detected 
radiation which was of a wavelength other than the characteristic Ka rays from the 
copper target. 

A calcite cleavage slip was mounted on the spectrometer and the x-ray spectrum 
was investigated, using as a detecting device in turn a Geiger counter and a 
proportional counter. The latter was used in conjunction with a high gain linear 
amplifier and -a_ single-channel pulse analyser adjusted so that only quanta 
corresponding to a narrow band of wavelengths centred on the Ka line were detected. 
The resulting curves are shown in fig. 1. Curve A is the intensity distribution 
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observed with a Geiger counter and curve B that with the proportional counter, 
a nickel filter being mounted in front of the counter in each case in order to reduce 
the intensity of the Kf line. Curve C is the intensity distribution at the x-ray tube 
window derived from curves A and B by correcting for the wavelength variation 
of the absorption in the intervening air, the nickel filter, the counter window and the 
counter gas. All three curves are also corrected for the variation of the reflection 
coefficient of calcite with wavelength (Compton and Allison 1935). For the sake of 
comparison the curves are plotted so as to make the intensity in the Ka line the same 
in each case. The total intensities in the two characteristic lines and in the continuous 
spectrum are compared in the table. It should be noted here that the sealed-off 


. Spectrum Geiger counter Prop. counter 
o 100 100 100 
B i353 0:09 0-067 
Continuous 18:5 3°5 0-9 


copper target x-ray tube was almost new and that it was operated at 30 kv d.c. In 
an older tube in which the target and the windows have become contaminated by 
tungsten, or in a self-rectifying tube operated from an a.c. supply, the ratio of energy 
in the continuous spectrum to that in the Ka line would be expected to be higher. 

The response of the conventional argon-filled Geiger counter to x-radiation 
of different wavelengths was calculated from published data (Internationale Tabellen, 
1935) on the absorption coefficients of the filling and of the window (curve A, fig. 2). 
It will be seen that the addition of a f-filter materially improves this response curve 
for the case of the copper lines (curve B). The curve, however, still shows appreciable 
sensitivity to the first harmonic, A/2 of CuKa, which may give rise to spurious effects. 
This sensitivity relative to the Ka wavelength cannot be appreciably reduced by the 
use of filters. A slight improvement can be effected by operating the x-ray tube 
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at a lower voltage with consequent loss of available intensity. 'The short wavelength 
limit of the continuous spectrum is marked in fig. 2 for two different excitation 
voltages. 

While the employment of filtered radiation combined with operation of the tube 
at a comparatively low voltage is adequate for many purposes, a very great 
improvement is achieved by substituting the proportional counter with appropriate 
circuits for the Geiger counter. Its response to any particular characteristic wave- 
length employed can be made as sharp as required. However, since the bias curve 
shows a spread about the centre of the selected channel, very sharp discrimination 
involves considerable loss of sensitivity to the Ka radiation. Fig. 2 curve C shows 
the response curve without nickel filter of the proportional counter used for a channel 
width of 3 volts when the centre of the channel at 28:5 volts corresponded to the 
Ke line, the external amplification being 10°. Under these conditions approximately 
50% of the Ka quanta absorbed in the counter were detected. This amount of 
discrimination reduces the sensitivity to the 4/2 harmonic to less than 5% of that 
for the Ka line itself. The non X-ray background of the counter is then about 
7 count/minute. 

The response curves for the counters may be compared with that for double- 
emulsion non-screen x-ray film shown in fig. 2 curve D, calculated on the basis of 
data given by Bromley and Herz (1950). 

The proportional counter used in these experiments was specially constructed 
for this purpose by Messrs. 20th Century Electronics Ltd. It had a length of 10 cm 
and a diameter of 5 cm. Its cathode was Aquadag coated on the glass walls. Radiation 
entered the counter normal to its axis through a berosilicate glass bubble window. 
The counter was filled with argon to a pressure of one atmosphere, about 8% of 
ethyl alcohol being added to stabilize the gas amplification. "The cathode was operated 
at —1750 volts, the 0-005-in. diameter tungsten wire being connected directly to the 
first grid of the pre-amplifier. The gas amplification under these conditions was about 
10°. The stability of the device when operated with a well-regulated power supply 
was in every way comparable with that of an externally quenched Geiger counter. 

In those cases where monochromatization by crystal reflection is still necessary, 
a proportional counter offers great advantages over a Geiger counter in that its 
sensitivity to the A/2 harmonic of the characteristic radiation can be made negligible 
and in that its background counting rate is about one-twentieth of that of a Geiger 
counter of comparable x-ray sensitivity. 


The Royal Institution, U. W. ARNDT. 
21 Albemarle Street, De Ra Rinnye 
London, W. 1. 


15th November 1951. 
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REVIEWS OF BOOKS 


Advanced Dynamics, by E. Howarp Smart. Vol. I, pp. xi+419; Vol. I, 
pp. xi+420. (London: Macmillan, 1951.) 40s. each volume. 


The two volumes of the work under review embody the substance of lectures given by 
the late Mr. E. Howard Smart, sometime Head of the Department of Mathematics at 
Birkbeck College, and they are obviously the work of a man of wide and varied experience as 
ateacher and examiner. One of the most useful features of the work is the range and number 
of examples which are solved in detail to illustrate general principles and theorems; these 
examples form a substantial proportion of the text, and they make the work very valuable to a 
student who is forced to grapple with the Miticuisies of dynamics without much help from a 
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teacher or supervisor. ‘The two volumes could, for example, be recommended without 
hesitation to a physicist who feels his lack of knowledge of the subject, and who is, at the 
same time, forced to remedy this state of affairs by his own unaided efforts. 

The scope of the work may be judged from the headings of the chapters. Volume I : 
Introductory, Fundamental dynamical principles, Rectilinear motion of a particle, Accelera-. 
tions parallel to rectangular axes, Radial and transversal accelerations (central forces). - 
Free motion under the law of inverse square, Constrained motion in two dimensions,. 
Motion in a resisting medium, Motion of a chain in two dimensions, Motion of a particle in. 
three dimensions; Volume II: Moments and products of inertia, Plane kinematics, 
D’Alembert’s principle and the equations of motion, Two-dimensional motion of a rigid 
body (finite forces), Two-dimensional motion of a rigid body (impulsive forces), Dimen- 
sions and dynamical similarity, Motion in three dimensions, Lagrange’s equations, Motion: 
under no forces, Motion of a top or gyrostat, Hamilton’s equations : general theorems on 
impulses, etc., Theory of small oscillations. 

The work is obviously written from the standpoint of a mathematician and mainly, one 
suspects, with an eye on the syllabus of certain specific examinations. For a physicist, one 
result of these restrictions is a weakness in the treatment of those parts of the subject where 
physics and applied mathematics meet, and by way of illustration one may quote the dis- 
appointing chapter in Volume II on Dimensions and dynamical similarity. A second result: 
is the omission of certain topics which are important to physicists; for example, whilst 
there is an exhaustive treatment of the small oscillations of a conservative system about a 
configuration of equilibrium, there is no mention of oscillations of finite amplitude or of 
non-linear oscillations. 

On opening these books one is immediately struck by the absence of Clarendon type, due- 
to the fact that the subject is treated almost entirely by resolution into Cartesian components, 
vector methods being ignored. This is not the place to marshal the case for and against: 
the vector treatment of dynamics; it is sufficient to record that this work offers neither 
apology nor explanation for omitting the use of vector methods, and, as a physicist, the: 
reviewer is not unsympathetic to this standpoint, if it were only because vector equations 
have to be translated into Cartesian language before they can be used in experimental. 
physics. R. M. D.. 


Integraltafel. Erster Teil: Unbestimmte Integral ; Zweiter Teil: Bestimmte 
Integrale, by W. GROBNER and N. HorreiTer. Pp. viii+166 ; vi+204. 
(Vienna and Innsbruck: Springer-Verlag, 1949 ; 1950.) 31s. ; 41s. 


These two volumes should find their way into most reference libraries. 'The authors. 
have taken advantage of earlier tables of integrals in making their compilation, and both the 
results taken from these and new results have all been independently checked. 

Volume I contains indefinite integrals, which are presented in three groups according 
to the form of the integrands. ‘These are rational functions, irrational algebraic functions. 
and transcendental functions. In appropriate cases alternative forms are given, for 
example in ascending and descending powers of the variable respectively. 

Volume IT contains definite integrals listed according to the same arrangement as that 
adopted for indefinite integrals in the first volume. 

The formulae and text have been photographed from a hand-written copy; although 
in German the script will be easily legible to the English reader. It is almost inevitable 
that some errors should occur in a collection of this size: the authors themselves quote 
three for volume I at the beginning of their second volume. 

Every specialist will note some omissions from his own field, but this is understandable 
and needs no excuse. Nevertheless it is a little surprising to note, for example, that no- 
indefinite integrals involving Bessel functions (of which there is a number of importance): 
are listed. 

These criticisms will not detract from the obvious worth of a collection of this kind. 


Both the authors and the publishers are to be congratulated on it. H. H. HOPKINS. 
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Dielectric Breakdown of Solids, by S. WuireHEaD. (Monographs on the 


Physics and Chemistry of Materials.) Pp. xv+271. (Oxford: Clarendon 
Press, £951.) 25s, 


About twenty years ago Dr. Whitehead published his first treatise on dielectric 
breakdown. Since that time so much progress has been made both experimentally, under 
the leadership of Von Hippel and of the author and his colleagues, and theoretically, mainly 
by Von Hippel and Fréhlich, that the author has now produced an entirely new book on 
the subject. 

Dielectric breakdown is a phenomenon of great practical importance and it has great 
theoretical interest. ‘There are at least three distinct types: intrinsic breakdown takes 
place almost instantaneously if the electric field reaches a level at which an electronic 
equilibrium is no longer possible; thermal breakdown is due to the heat developed in the 
dielectric as Joule heat or by dielectric losses; and discharges outside the dielectric or in 
small cavities inside can weaken it chemically or thermally to the point of breakdown. 
It is the burden of Dr. Whitehead’s book to set out these types, to show and to illustrate 
how the dielectric strength in its different aspects is measured, to summarize the results 
obtained and to explain the theoretical background and the detailed theories and finally 
to assess the importance of the different aspects of breakdown and the future trends of 
practice and research. In addition a discussion is given of the variation with time of 
dielectric strength due to diffusion and slow chemical processes. 

This task is even more formidable than it may appear at first sight. ‘The author tries to 
reduce to a compass of about 250 pages the contents of about 200 original papers, the 
scope of which includes subtle experimenting and fundamental theoretical analysis of 
electronic processes in dielectrics; and as original papers are, these days, usually already 
condensed almost beyond the limit of comprehensibility, the book makes very heavy 
demands on the reader’s acumen. Details concerning experimental arrangements are 
scanty and the reader needs a great deal of imagination to get an idea how measurements 
are actually made and how the results given in diagrams are obtained. Even the plates 
showing photographs of faulty dielectrics are barely explained. ‘The mathematical theories 
are given an appreciable portion of the space available; but the story is so involved that 
only readers with a good knowledge of the original literature may be able to follow. 

The book will be warmly welcomed by those specially interested in dielectric 
breakdown as a clearly laid out summary of the subject written from a very high level. 
It is further an important book of reference for scientists interested in electronic processes 
in solids. It will give the reader readily an idea of the state of the subject, to the extent 
to which a much shorter review article could. But it will yield its full contents only to 
those who are prepared for a great deal of additional study. W. E. 


Tensor Analysis for Physicists, by J. A. SCHOUTEN. Pp. x+275. (Oxford: 
Clarendon Press, 1951.) 30s. 


This is no easy introduction to the subject of tensor analysis. The outlook of the 
author is that of a mathematician, and, though he is trying to write for physicists, he cannot, 
in the nature of things, share their outlook. As a consequence, the first part of the book 
will appear to many to be forbiddingly difficult, though in fact it only wants perseverance 
to grasp any of the ideas used. The later part of the book, where the general methods of 
tensor analysis are applied in turn to elasticity, classical dynamics, relativity and quantum 
mechanics, will be found in some ways easier to read than the earlier parts. 

Much of the difficulty felt in reading work expressed in tensorial form is due to the 
unfamiliarity of the notation; Professor Schouten uses very clear and unambiguous 
symbolism, but one which makes considerable demands on the attention of the reader, 
who has to remember, for example, to distinguish between roman and italic numerals, 
as well as letters. 

The book is exceptionally well produced, and is a valuable reference work to be included 
in any serious library of physical works. Ifo 186. EMA 
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ABSTRACTS FOR SECTION B 


Water Bells, by F. L. Hopwoop. 


ABSTRACT. The author shows that when water is discharged through an annular slit 
above a plane water surface, the sealed “ water bells’ thus formed exhibit some novel and 
remarkable properties not hitherto described. The primitive water bell is a dome-shaped 
bubble whose dimensions increase with increased rates of flow. If this is perforated with a 
finger the maximum diameter of the bell suddenly doubles itself. On continuously reducing 
the flow the expanded bubble contracts and assumes an alternating sequence of stable and 
semi-stable forms of great beauty. All these, when perforated, contract slightly. The 
semi-stable forms have the general appearance of a hyperboloid surmounted by a saucer- 
shaped depression with the annular slit at the bottom. On inflating an expanded bubble by 
a slow stream of small air bubbles, semi-stable forms similar to the above are produced, 
but possessing an additional inflection in their contours. During these changes the 
maximum differences of pressure do not exceed one-tenth of a millimetre of water above or 
below atmospheric pressure. ‘The author shows how the apparatus can easily be modified 
to produce enclosing water bells. It may also be adapted to produce sonic and ultrasonic 
underwater vibrations. 


The Use of Current Sheets in the Design of Magnets to give Bounded Fields of 
Required Form, free from Edge Distortion, by H. O. W. RicHarpson. 
ABSTRACT. The use of current sheets to produce static magnetic fields of required 
form, with freedom from edge or fringing distortion and with no stray flux, is discussed. 
Two types of magnet winding are considered. The first, or field- -terminating type, encloses 
a tube of force and serves as a sharp boundary to the field. The second, or iron-shrouded 
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absorption (medium Hilger, zenith plate). 
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Fig. 1. Pattern from a KBr cleavage face, upon which a thin layer of silver chloride has been 
evaporated. KBr [001] azimuth. 


lig. 3. Oriented TIC] on a KBr cleavage Fig. 4. TICI layer on a mica cleavage face. 
face. 'T'IC] has an abnormal NaCl type Mica [010] azimuth. 


structure. KBr [001] azimuth. 
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type, is wound against the inner surface of an iron shroud and, with a suitable distribution of 
turns, can generate any required shape of field in the iron enclosure. In both types stray 
field can be avoided if the iron return-path is infinitely permeable, or if it is covered with an 
“anti-reluctance’ winding. 

A discussion is given of (a) the uniform field, (b) the field whose lines of force are coaxial 


circles, (c) the synchrotron field, and (d) the prolate spheroidal field, whose lines are confocal 
ellipses. 


The Extension of the Electrolytic Tank Method to the Study of Magnetic Fields due 
to Iron-clad Current Sheets in Three Dimensions, by H. O. W. RICHARDSON. 


ABSTRACT. An application of the electrolytic tank method to plot the magnetic fields of 
a current sheet is described. In general, three-dimensional fields due to currents can be 
plotted, provided that the current is close to the surface of a piece of permeable iron. 


On the Torsional Vibrations of Long Chain Molecules, by B. Sz1GrTI. 


ABSTRACT. Simplified models are discussed for long chain molecules in which all 
units are equal except for one which is heavier than the others. Certain general conclusions. 
are reached for the vibrations of such molecules. The results will be applied to the out-of- 
plane or torsional vibrations of long chain ketones in a following paper. 


Fluorescence of Dyestuff Solutions: Viscosity and Temperature Effects, by 
N. R. Tawpbe and N. RAMANATHAN. 


ABSTRACT. The effect of viscosity and temperature of the medium on the polarization 
of the fluorescent light emitted by some dyestuff solutions when excited by linearly 
polarized monochromatic light is investigated by the visual and photoelectric methods. 
The variations in the polarization caused by temperature are found to depend on the 
concentration and molecular weight of the dyestuff. 


Optical Transmission of Certain Colloidal Solutions under Relaxation Conditions, 
by G. KLEIN. 


ABSTRACT. Formulae are derived for the variation in time of transmission of light 
through dilute colloidal solutions of uniaxial particles which are initially orientated by an 
external electric or magnetic field from the instant the field is removed. ‘'I'wo cases are 
considered : (i) particles opaque and large compared with the wavelength of light, and 
(ii) metallic particles small compared with the wavelength. 


An Interferometric Method of Studying Local Variations in the Refractive Indices 
of a Solid, by R. C. Faust. 


ABSTRACT. Many small transparent solids possess local variations in their refractive 
indices that cannot be detected by the standard immersion methods. This difficulty has 
been overcome by placing the immersed solid between the two parallel, metal-coated 
plates of an interferometer illuminated at normal incidence by a parallel beam of white 
light. The immersion liquid is so chosen that its dispersion curve intersects that of the 
solid at some wavelength within the visible spectrum. An image of the immersed specimen 
is then projected on to the slit of a spectrometer, whereupon one can observe coloured 
fringes of equal chromatic order (channelled spectra). ‘These fringes generally suffer 
a displacement on crossing from the liquid to the liquid-solid region of the interferometer ; 
this displacement does, however, vanish at that wavelength at which the refractive index 
of the liquid equals that of the solid. The method is equally applicable to both isotropic 
and birefringent specimens, the index measurements under favourable conditions being 
accurate to within +0-0001. 

The paper is illustrated by several examples, particular reference being made to the 
possibility of studying the skin effect in rayon fibres. 


80 Abstracts for Section B 


Optical Properties of Tellurium in the Infra-Red, by 'T. S. Moss. 
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ABSTRACT. The refractive index and absorption constant of tellurium layers have been x | 


measured in the near infra-red. The refractive index has been measured over the wave- 


length range 3 4 to 11. The value extrapolated to long wavelengths is found to be 4:8; 


this is much higher than any other known refractive index. Absorption measurements were — 


made for wavelengths 1. to 3 yz: the absorption constant is 2-2 x 10° cm~? at 1:0 p, falling 
rapidly with increasing wavelength. Transmission measurements on bulk specimens are 
possible at wavelengths greater than 3:3 p. 


The Diffraction of X-Rays by Distorted-Crystal Aggregates: III—Remarks on the 


Interpretation of the Fourier Coefficients, by J. N. EasTaBrooxk and A. J. C. 
WILSON. 


ABSTRACT. - The line profile of the Ath order reflection from a set of lattice planes is 
known to be given by a Fourier series with coefficients A(h, m)=N(m)J(h, m). The 
function N(m) depends on the particle shapes and sizes, and J(h, m) depends on the state 
of strain. Neither can be completely determined from the experimental Fourier coeffi- 
cients for a single value of h. They have, however, different behaviours for small m, and 
A(h, m) has an unambiguous interpretation in this region only : its initial slope gives the 
mean particle thickness, and its initial curvature sets a lower limit to the mean-square strain. 
The practical value of this interpretation is subject to two limitations : (i) the initial shape 
of the A(h, m) curve depends on the ‘ tail’ of the line profile, (ii) if the crystals are large 
enough to contain several regions of compression and extension, probability theory shows 
that J(m) is likely to be approximately proportional to exp{—const. |m|} for m large. 
Extrapolation to small m could then tend to give fictitiously small values for both particle 
size and strain. 

The problem of obtaining the greatest possible amount of information from observed 
line profiles of all orders is discussed. “Usually only a few orders can be measured, and the 
information obtainable is less, but in three special cases the Fourier coefficients oft two lines 
can yield all information obtainable from the full set. 


A New Determination of the Luminance Factor of Magnesium Oxide, by J. S. 
PRESTON and G. W. GorDon-SMITH. 


ABSTRACT. A simple, direct, and absolute method of reflectometry is described, 
avoiding the use of the inverse square law, and overcoming several of the drawbacks of the 
selenium photovoltaic cell which is used for the measurements. These advantages are 
secured through the use of calibrated apertures and certain integrating properties of the 
cell. The result obtained for the luminance factor for white light of an MgO smoked 
layer 1 mm thick was 1:01 for normal incidence and 45° view. 


